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ABSTRACT 


Numerical methods to determine the flow field around blunt bodies 
traveling at hypersonic speeds in planetary atmospheres are presented. 

The state and motion variables and the convective and radiative heat 
transfer are determined for both equilibrium and non-equilibrium flow 
models. The flow model is chosen to be applicable to conditions en- 
countered in typical out-of-orbit and hyperbolic Martian entry tra- 
jectories where maximum stagnation point convective and radiative heat 
transfer occurs. For the conditions of interest, the coupling between 
the inviscid region, the viscous region, and the radiative transfer is 
significant, and is included through simplified engineering approximations 
Numerical results are presented for equilibrium and non-equilibrium flow 
over a sphere-cone having a base diameter of 12 ft, and moving at 19,600 

ft/sec in a 70% N„, 30% CO atmosphere at an ambient density of 10 ^ 

3 LI 

slugs/ft . Results are also given for non-equilibrium flow over a sphere- 
cap at the same flight condition. 
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SECTION L 
INTRODUCTION 


The flight conditions and the planned vehicle shapes for Martian entry 
pose a number of problems m making satisfactory flow field predictions. 
Because of the low ballistic coefficient vehicles, peak heating occurs at 
low ambient density, which leads to a relatively thick viscous region. 

The high emissivity of CO^'N^ mixtures leads to significant radiative 
heat transfer. Because of the low ambient density, much of the gas in 
the shock layer is out of equilibrium and produces larger radiation emission 
than would be produced if the gas were at equilibrium. Finally, because of 
the low density and high emissivity, the radiation losses from the shock 
layer gas lead to a reduction m gas enthalpy. To carry out the desired 
flow field calculations, the flow model must include all of these physical 
effects, but it should not become unnecessarily complicated. Thus we have 
chosen the simplest methods of calculation that retain the essential 
physical information. 

A survey of applicable theoretical and experimental low Reynolds number 
shock layer studies indicates a wide variation m the method of solution, 
the gas models, and the assumed free stream conditions with a corresponding 
large variation m the predicted magnitude of "rarefied regime" effects. 
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At present one of the most general, and yet definitive studies of rarefied 
gas regimes is that given by Goldberg.^ In this study, non-similar hyper- 
sonic continuum flow field sphere-cap solutions were obtained for a flow 

, 2 
model which is applicable through the incipient merged layer regime. An 

examination of the results of References 1 and. 3 indicates that for the 
shock layer Reynolds numbers of interest m this study, viscous and conduc- 
tion effects due to velocity and temperature gradients at the shock and 
£n the inviscid shock layer can be neglected. Thus, Rankine-Hugoniot 
shock relations and inviscid shock layer equations can be used to calculate 
the flow field. It also is assumed that the vehicle wall is- cold m this 
study and, thus, the boundary layer displacement thickness is small and the 
inviscid shock shape and body pressure distribution are not affected 
appreciably by the boundary layer. 

For the flight conditions given m this study, the low Reynolds number 
1 3 

results 5 predict greater convective heat transfer rates than are obtained 
from the classical discontinuous shock - inviscid shock layer - thin 
boundary layer theory. This increase in heat transfer is primarily due 
to the shock layer vorticity induced by the curved bow shock. The effect 
of vorticity can be incorporated in conventional boundary layer theory, 
however, as long as the proper boundary layer edge velocity conditions 

.1,3 

are used. ’ 

In computing radiative heat transfer it is apparent from shock tube data 
that for the flight conditions under study much of the gas m the shock 
layer is out of equilibrium. However techniques for computing non- 
equilibrium radiative heat transfer are not yet well developed. Thus it 
is desirable to carry out equilibrium calculations as well; although an 
equilibrium flow model is less realistic here, the accuracy of the 
numerical results is greater. In Section 2 below, the equilibrium flow 
model is described and in Section 3 the non-equilibrium model is set 
forth. The numerical results are given in Section 4. 
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SECTION 2 
EQUILIBRIUM FLOW 


2.1 INVISCID ANALYSIS 

It was originally planned to use the first approximation method of integral 
relations computer program to determine the mviscid flow properties. A 
number of modifications in the analysis were needed, however. As these 
were investigated, it became apparent that several difficulties could not 
be overcome satisfactorily. The major problem with the first approximation 
integral relation solution is that, at the low specific heat ratios which 
are representative of the equilibrium case, the 60° cone shock layer thick- 
ness is significantly overestimated. This difficulty is illustrated further 
in Section 4.1. It was also found that the stagnation point velocity 

4 

gradient was not determined accurately by the one strip integral method 
. and that attempts to infer this quantity from the numerical results were 
unsuccessful. Furthermore integration of the inviscid flow equations 
through the transonic region was found to be very difficult and could not 
be carried out satisfactorily. Thus an alternative and, fortunately, more 
accurate method was applied to this problem. This method was under develop- 
ment at the time the present program was planned. 
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The method has been developed here for the computation of either planar 
or axisymmetric blunt body flow fields. For steady flow fields the formula- 
tion of the problem is direct in the sense that the body is given and the 
shock wave computed. Steady flows are computed as the large time limits 
of transient flows. The time dependent equations of motion are solved by 
a finite difference technique. The stable, explicit finite difference 
approximations are determined by a method proposed by Godunov, Zabrodin 
and Prokopov.^ 

An initial condition in the form of velocity components, density, pressure 
and internal energy plus kinetic energy is specified on a mesh covering , 
the region of integration. This initial data is m theory arbitrary, but 
in practice is specified by a steady, constant y flow field. The amount 
by which the data does not satisfy the steady equations of motion determines 
the time derivative of the flow field. The time derivative then determines 
new values of the flow variables in each cell of the mesh. In this manner, 
the flow field evolves in time, subject to the conservation laws, from the 
specified initial conditions. 

The boundaries of the mesh are dependent on the method by which the bow 
shock is treated. It has been common practice m finite difference 
computations to ignore the existence of shock waves and leave their treat- 
ment to the finite difference scheme itself. The computed shocks are 
represented by smooth but rapid transitions of flow variables over a 
linear distance of nominally 5 to 10 mesh points. To adequately model 
the flow variations, 5 mesh points should constitute a small portion of 
the shock layer width, and commonly 50 to 100 points would be required 
between the body and the upstream mesh boundary at which a uniform stream 

condition is applied. Smeared shock calculations of this type have proved 

£ 

to be quite feasible for blunt body flows, Bohachevsky has used Lax's 
scheme,^ Burstem^ used the Lax-Wendroff scheme,^ and Emery^ compared 
Rusanov's scheme’*"^ with those of Lax and Lax-Wendroff. 
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The alternative employed in this code is to retain the bow shock as a 
discontinuous boundary of the region of integration. To achieve this 
distribution of computational points, it is necessary to introduce a mesh 
that always exactly fills the shock layer. The cells in the mesh expand 
and contract as the standoff distance between the shock and body increases 
or decreases as a function of time and position. The upstream boundary of 
the mesh maintains the full jump in flow variables from the free stream to 
the shock layer values immediately behind the shock. This type of mesh has 

been widely used in the blunt body problem; m the time dependent form, see 

2 5 12 

Hayes and Probstein, Godunov, Zabrodin and Prokopov, Moretti and Abbett. 

In addition ' to the bow shock wave, the boundary of the region of integration 
is composed of the x-axis, the trace of the body m the xy plane, and a 
surface entirely embedded in the supersonic portion of the flow, across 
which there is no upstream influence. 

Given the upstream and body portions of the mesh boundary, the internal 
mesh is constructed in the following manner. Straight line rays are 
selected which intersect both the body and the upstream shock boundaries. 

The ith ray makes an angle 9^ with respect to the x-axis and intersects 
the axis at the position x = 5 as illustrated in Figure 1. The rays 
divide the shock layer into discrete intervals along the body. The first 
of these rays is coincident with the x-axis, and the last forms the down- 
stream boundary, immersed m supersonic flow. 

The interval on each ray between the interaction with the upstream and 
body boundaries ls'divided into a fixed number of submLervals , say J. 

Let . represent the distance from the x-axis measured along the ray 
to the intersection of the ray and upstream-shock boundary, and ^ 
denote the similarly measured body-ray intersections. Then the lengths 
•b. , and b . denote the inner and outer end points of the ith subinterval 

j+l,i j,i 

in from the bow shock boundary. 
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FIGURE 1. MESH GEOMETRY FOR SHOCK LAYER AND 
CELL j+1/2 , l+I/2 
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When the points -L. for fixed j are connected by sLraight line segments, 
a polygonal interpolation curve running along the shock layer is obtained. 
These interpolating polygons and the rays divide the shock layer into a 
mesh of quadrilateral cells upon which the difference equations are 
constructed. 


The cell bounded by the ith and i+lst rays and the jth and j+lst transverse 
interpolating polygon is denoted by j+%, i+%. Similarly, the interval on 
the ith ray bounded by the intersections with the jth and j+lst inter- 
polating polygon is denoted by j+%,i. Finally the segment of the j th 
interpolating polygon, j, i+% is that segment running between the ith 
and i+lst rays, see Figure 1. 


2.1.1 FINITE DIFFERENCE EQUATIONS 

The differential conservation laws will be used in the form 


f + P + Q +y _1 R=0 
t X y 


( 1 ) 


The quantities f, P, Q, R are the following vector valued functions 


f = 


P 

pu 

Pv 

-_ pe u 


( 2 ) 


P(f) = 


pu 

2 , 

pu +p 


puv 

pu£+up 


(3) 



Q(f) 


pv 

puv 

2 

pv -f-p 
pve+vp 


R(f) = k 


pv 

puv 

2 

pv 

pv(e + 


( 4 ) 


( 5 ) 


where k = 0,1 for plane and axisymmetnc problems respectively and 

2 2 

e = e(p,P) + %(u +v ) is the internal plus kinetic energy. 


The difference equations are generated by integrating Equation (1) over a 

cell 0(t), whose boundaries are time dependent due to the floating mesh, 

and from initial time t to t +', The integrated form reads 

o o 

t +T 
o 

fj f(x,y,t +T)dxdy = [T f(x,y,t Jdxdy -b f dt if f(x,y,t)V ds 

n(t a ( t > t aa(t) n 

o o o 


- J) (Pdy - Qdx) - fj ” R(x,y , t) dydx} (6) 

an(t) Q(t) y 

where BQ is the boundary of 0 and represents the component of the 

boundary velocity m the outward normal direction. This equation is to be 

evaluated to first order in smallness of time step r and altitudes h ,h 

X Y 

of the cell Q. The derivation of the finite difference equations is 
completed in Appendix A. 
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2.1.2 SOLUTION OF DIFFERENCE EQUATIONS AND EQUATION OF STATE 
An equation of state must be added to finite difference Equations (A13) 
to (A16) before a computation can be carried through. For equilibrium 

runs, an Ames polynomial fit to the equation of state of a 307<> C0„-707 o 

14 1 

^ mixture was employed. The polynomials provide the internal energy, 

e, and sound speed, a, as a function of p and p. However p and e are the 
thermodynamic variables determined at the new time by the conservation 
laws. Thus some iteration process is necessary to determine p accurately 
at each time. This iteration can be done simultaneously with the flow 
field's approach to steadiness in such a manner as to require no addi- 
tional evaluations of the polynomials. 


The initial thermodynamic data in the cells is specified by and P^°^. 

From this data the internal energy is determined from the Ames tables. 

Jo) 


The ratio r^°^ is defined by 


» . 


p (o)~(o) 


(7) 


At the nth time step, the quantities and r*' 11 ^ are 

s t 

determined from these same quantities at the n-1 time step m the 

following manner. The conservation laws will provide the quantities p 

, (n) _ (n-1) — (n-1) , (n-1) , , , 

and e rrom p , e and p . The new pressure will be given by 


(n) 


(n) (n-1) (n) (n) 

p v = r p e 

The Ames tables will give e^ n ^ as a function of p^ n ^ and p^ n ^ and the 


( 8 ) 


ratio r 


(n) 


follows as 


.(n) = 


» 


p (n)-(n) 


(9) 


When steadiness of the flow field is achieved p*' 0 ^ P^° ^ 

, -(n-1) -*(n) „ „ . — (n) (n) , , 

and e = e From Equation (o) e = e and hence a consistent 

answer is achieved on convergence. The convergence of this scheme depends 

upon the polynomial coefficients in the Ames tables. Qualitatively 
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convergence is achieved since the quantity r varies very slowly with any of 
the thermodynamic quantities. Over the entire flow field r will vary by 
about 207» and by no more than 17. over a time step. Thus the p^ n ^ 
predicted by Equation (8) will be quite accurate for use in the conserva- 
tion laws and e^ n ^ and e^ n "* will not differ appreciably. 

2.1.3 FLUXES ACROSS CELL BOUNDARIES 

The flow quantities on the cell boundary segments internal to the mesh 

are given in Godunov's scheme by the solution of a related one dimensional 

unsteady flow problem. The one dimensional problem is constructed using 

the pressure, density and normal component of the fluid velocity relative 

to the cell boundary velocity in the two cells adjacent to the segment. 

These values are used to define the initial uniform states in a Riemann 

problem, the general discontinuity separating two uniform fluid states. 

This discontinuity evolves m time through left and right running 

disturbances separated by a contact discontinuity. The full solution for 

13 

the properties is prescribed in Courant and Friedrichs. 

As numerically demonstrated by Godunov, stability is maintained, even 
when computing through a shock, when the solution of the Riemann problem 
is linearized with respect to small differences across the discontinuity. 
Then any rarefaction disturbances have negligible width and the cell 
boundary lies m one of four constant state regions: (1) to the left of 

the left running disturbance, (2) between the left running disturbance 
and the contact surface, (3) between the contact surface and the right 
running disturbance, (4) to the right of the right running disturbance, 
see Figure 2. Note that the t axis is the position of the cell boundary. 

If the disturbances propagate so that the t axis is located in regions (1) 
or (4), the cell boundary lies in undisturbed fluid. The pressure, density 
and velocity components on the boundary are equal to the quantities in the 
cells corresponding to the initial left or right states. If the t axis 
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CONTACT 

DISCONTINUITY 


/ 


li 

LO 

'w' 

(3) / 


n,c 

RIGHT 

= p (3) 

yS 

RUNNING 

(3) 

yS 

FRONT 


^ "i* 

w 

n,r 

p+ 

P+ (4) 



IN FLUID STATE 



is located in regions (2) or (3), the cell boundary lies behind a 

disturbance, and the pressure and normal velocity component are determined 

by expansion of the solution in Reference 13. Denoting the average of 

the initial right and left state by bars and the difference of the left 

from the right state by A's, permits the pressure and normal relative 

velocity components, w , in regions (2) and (3) to be written as 

n * x 


? (2) = p(3) = p . ^ /2 +• 0(u 2 ) 


= w = w - £ip/2m b 0(A 2 ) 

n,r n,r n,r 


w 


( 10 ) 

(ID 


m = p a b 0(A ) 


( 12 ) 


where a is the sound speed evaluated at p, p. The entropy s is constant 
across the disturbances 


s 

s 


( 1 ) = 
(3) = 



(13) 

(14) 


( 2 ) 

The pressure on the cell boundary has the value p , and the normal 

component of the fluid velocity relative to the boundary velocity is 

( 2 ) 

w . The density on the cell boundary is determined by the known 

n , r 

values of pressure and entropy, and the fluid velocity component tan- 
gential to the cell boundary has a value equal to those in regions (1) 
and (4) depending on whether the t-axis lies to the left or right of the 
contact discontinuity. 


2.1.4 BOW SHOCK BOUNDARY CONDITIONS 

At the beginning of a time step, say t Q , each bow shock segment separates 

a cell containing a uniform state of fluid from the free stream. Let the 

+ + + + 

state in the cell bordering the bow shock be denoted by p , p , w , w , 
b b n t 

where w , w are the components of the fluid velocity resolved normal 

and tangential to the bow shock segment. The jump in the properties p + , 

p , w to the free stream values p m , P . q can be viewed as a Riemann 

* 1*2 1 co J co 1 1 co n 
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problem in a coordinate n normal to the segment of the bow shock, with 
n > 0 directed downstream and n < 0 upstream. The left running disturbance 
must be identified with the actual bow wave. The right running dis- 
turbance carries the information of the bow wave's motion during the 
ensuing time interval into the first cell. 


For the bow shock, if 


= ^P cd -P~)( 1 /P _ - 1 -/P~ cd )3 2 


then the momentum equation is 


w m + p = w ,m 
n ncd 


cd 


and the Hugoniot relation is 


e cd " e = ^ (P cd + P ) (x/p ' L/p cd ) 


(15) 


(16) 


(17) 


The right running disturbance on the other hand is weak and travels with 
the sound speed, a + , into the interior of the cell. Then 

m + = p + a + (18) 


is known explicitly. The momentum balance across this disturbance is 


h , 

w ,m + 
n , cd 


cd 


h + + 

w m + p 
n 


(19) 


Equations (15)-(19) can be iterated to determine the pressure and normal 
velocity on the contact surface, p , and w , which determine in turn 
the quantities R, U, V, P, E to be used in Equations (A13) through (A16) . 

If the cell under consideration is 3/2, x+h the fluxes are to be determined 
from 
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R 


l,i+% 


cd 


^n"* 1 j i-b'j ''cd 

<Vi,« ' q - COOT i,iH 


V+* ■ p 


cd 


( 20 ) 


where o t is the angle of the bow shock segment 1, i+k relative to the 
I , i+% 

free stream direction. The velocity of the bow wave, V , can be found, 
after m is known, from the mass conservation principle 


P«,(<ho n - V ) = m 


(21) 


2.1.5 MESH DEFORMATION DURING TIME STEP 

Since the upstream boundary of the mesh is the only portion of the boundary 
that moves, the deformation of the mesh during the time increment follows 
from the analysis of the preceding section. If is the velocity of the 
intersection of the shock with the ith ray inwardly directed, then the 
mesh deformation will be determined by the differential equation 

d-C 

— -1 - V. (22) 

After each segment of the bow shock is allowed to move parallel to itself 
during the increment T, the shock, though continuous initially, would not 
end up continuous. What was neglected above was the angular velocity of 
rotation of the segment, d .. Let V ., denote the shock velocity for 
the shock segments just above and below the ray i. The V x of course 

l+'2 

are the V determined previously for the two cells. From Figure 3 it is 
seen that can be expressed m two ways one relative to the lower ray, 
the other relative to the upper. It follows that 

v _ V i+h " ”l,iH--2 d l>i+% 

1 " Sin(CT l,i+^ 6 i ) 2 sin(CT 1 ,i ± % +e i ) 
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t-ofM 


where the lower or upper signs are used consistently throughout the 
expression and d^ are the lengths of the bow wave segments above and 
below ray i. If the expression with the upper sign is multiplied by 




and the expression with lower sign by 


sin (a, i+9.) 


“l,i-f% -% 

and the resulting expressions added, then V. can be expressed as 


V. f 

i * 


V. -d. , , + V ; d. . , 

i+% l,i+% i-% l,i--g 


i ' d, . , sin(a, , ,+9.) + d, . ,, sin(o, ,+9.) 
l,i-h l,i+iz l l,i+% l,i-% i 




“‘Vd 1 


(23) 


where only the order of the coefficient of (d ,-d. x ) is indicated. 

I j 1 +^ L ,x--i 

The first term on the right is seen to be 0(1) while the last term is 

2 2 

0(d ,) = 0(h ) and is to be neglected in a first order scheme. Thus 

1 J 1+T2 

the motion of the shock node on each ray is determined by the motion of 
the segments adjacent to that ray. The shock shape after T is found by 
connecting the new nodal locations. The internal mesh deformation is 
subsequently found by interpolation on each ray with respect to the new 
standoff distance found from Equations (22) and (23). 


2.1.6 BODY BOUNDARY CONDITION 

The cells J+%, i+% bordering the body satisfy the condition of zero 

normal velocity on their body boundary J=l, i+%. Let the cell state be 
+ + 4 - + + 4 - 

denoted by p , p , w , w where w , w are the fluid velocity com- 

n t n t 

ponents normal and tangential to the body. Since the body boundary is 
stationary, the normal velocity condition is satisfied if a reflecting 

cell located immediately inside the boundary is imagined with state p = 

+ - + - + , + 

P , P = P , w t = w and w^ = w^ . 

+ 

Since m = -m by Equation (25), Equations (26) and (29) require that 

w , = 0 which is the zero normal velocity condition. The pressure at 
n,cd 

the wall boundary is 
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( 24 ) 


+■ + + 
p , = w m + p 
r cd n 

This Riemann problem features a contact discontinuity which remains 
coincident with the body boundary. 

2.1.7 EVALUATION OF THE ACCURACY OF THE INVISCID ANALYSIS 
The accuracy of the method given here to compute mviscid flows has been 
evaluated by carrying out several test calculations. In Figure 4, the 
surface pressure on a sphere in a M = 4.0, Y = 1.4 free stream is compared 
with two other numerical calculations. The results of the present code 
contain some scatter, but the deviation from Belotserkovski 1 s third 
approximation, scheme I, integral solution is less than 3%. A further 
confirmation of this estimate of the accuracy is that in the fine and coarse 
mesh comparisons of the blunt cone in an ideal gas the surface pressure on 
the cone differs by 27». 

A shock shape comparison with data is shown m Figure 5. On the conical 
portion the computed shock shape lies outside the data by less than 5%. 
However m the stagnation region it is difficult to make an accurate 
assessment from the photograph. In this region, comparison between the 
fine and coarse mesh computation of the blunt cone flow field again provides 
the error estimate. The coarser mesh is larger by 3.1% than the fine mesh 
value for y ~ 1.1. The density along the entire stagnation streamline 
for the coarse mesh lies about 5 % below the values for the f me mesh. If 
we assume a linear truncation error which is of order 3 to 57o when the mesh 
width is reduced by then the fine mesh results lie within 3 to 57 > of 
the actual values with no truncation error and the coarse mesh values lie 
within 67» to 10% of the actual values. 


The accuracy of the stagnation point velocity gradient was estimated during 


the checkout phase of the finite difference program by comparison with 

„ , . 15,16,17 . , 18,19 , . . . „ 20,21 

existing data, correlations, and numerical solutions. 

Spheres in an ideal gas stream with y = 1.4 at Mach numbers of 4.0, 5.0, 
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SURFACE TO STAGNATION POINT PRESSURE (P,/P 
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Time Dependent 
Solution 



FIGURE 5. BLUNTED CONE SHOCK SHAPE COMPARISON WITH SCHLIER5N. 

Moo= 9.0, 7 == 1 4 SCULIEREN FROM DATA OF NICHOLS AND 
NIERENGRATEN 82 



8.0, and 10.0 were computed as the check cases. For these runs, the shock 
layer was covered by a mesh with 7 cells across the layer and 6° increments 
along the surface. In Figure 6, a compilation of these results is pre- 
sented, in which it is seen that the stagnation point velocity gradient 
computed by the Aeronutronic program lies well within the scatter of data 
for Mach numbers less than 5.0. There is a serious scatter between data 
correlation and theory at Mach numbers greater than 5.0. However, the 
Li- Geiger correlation formula, the two Belotserkovskii integral relations 
solutions (second approximation, scheme I and fifth approximation, 
scheme III), Manglers inverse method calculation and the Aeronutronic 
Godunov type program scatter less than 57. at these larger Mach numbers. 

The Belotserkovskii and Mangier results are inherently the mast accurate, 
in which case, the Godunov finite difference program produces results about 
2-4% too low. 


For the 60° blunt cone shape, an ideal gas computation was performed to 
give initial conditions for the required equilibrium solution. Two meshes 
were used both of which had 7.5° intervals on the spherical portion while 
3 and 6 cells were taken across the layer. The free stream conditions 
were specified by a Mach number of 20.0 and y = 1.1. The coarse mesh 
computation gave 


iO 0.3 ■ - 366 


» V ■ 1.1 


while the fine mesh yielded 


— 0^), , = .372 
q ro os 0,6 


, Y = 1.1 


These results agree closely with the value predicted by the constant 
density solution for a sphere 


. (8 

qJ'os'OjlO v 3 P 2 J 


.365 
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STAGNATION POINT VELOCITY GRADIENT 


O KOROBKIN AND GRUENKWALD 15 
+ KENDALL 16 
A BOISON AND CURTIS 17 

@ EELOTSERKOVSKI 21 SCHEME I, SECOND APPROX 

SCHEME III, FIFTH APPROX 
90 

Xmangler 

LI AND GEIGER 19 

A SPHERES BY PRESENT PROGRAM 



FIGURE 6. STAGNATION POINT VELOCITY GRADIENT ON SPHERES 
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This agreement with the asymptotic sphere values reflects the weak 

upstream influence of the conical surface in the y - 1.1 case. Since a 
. o 

60 cone possesses an attached hypersonic shock layer solution, the only 

upstream influence of the conical portion is due to a narrow entropy layer 

at the surface. For a larger effective y or larger cone angle, this 

phenomenon would not occur and a serious divergence between the spherical 

20 

prediction and the actual computed value would arise. The Li-Geiger 
19 

correlation is in error at the large density ratios encountered here 
as it has the wrong a'symptotic limit as y - 4 1. Its predicted value is 
.315. 


The equilibrium computation provided a stagnation point velocity gradient 
of 


equil ’ ' 388 F 


= 6340 sec 


-1 


where the mesh is the coarse mesh in the preceding paragraph. This value 
corresponds to a constant density shock layer on a sphere with effective 
V = 1.125, in agreement with the ratio h/e computed in the flow field. 

The error in this computed value of the stagnation point velocity gradient 
is estimated at less than 5%. This error is based on the average 3 °L error 
between the y = 1.4 computations and those of Mangier and Belotserkovskii 
and the 1-2% difference noted between the fine and coarse mesh values. 


2.2 CONVECTIVE HEAT TRANSFER ANALYSIS 

The convective heat transfer calculations are complicated by the presence 
of low Reynolds number effects, significant quantities of C0^ m the 
boundary layer, and dissociation and ionization effects. The following 
paragraphs outline the analytical methods, including the above considera- 
tions, which were utilized in this study. 

In the discussion of the flow model in the introduction it was concluded 
that low Reynolds number effects are primarily manifested as vorticity 
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corrections to conventional boundary layer theory. At the stagnation 
point and on the spherical portion of the sphere-cone and the sphere-cap, 
the vorticity interaction effects have been calculated by using the method 
presented in Reference 3 and by utilizing the low Reynolds number sphere- 
cap solutions given in Reference 1. The method of Reference 3 is a solu- 
tion to the boundary layer equations in which the slope of the boundary 
layer velocity profile is matched with the slope of the inviscid velocity 

profile for equal mass flow rates. The input parameters required to ob~ 

2 

tain a solution are the body Reynolds number, R g = (du/ox)/|i o 

o 

and the normal velocity gradient, d(u/u Q ) /S(y/R^), These parameters are 
given by the inviscid flow field solution. It was found that the method 
of Reference 3 and that of Reference 1 both give essentially the same 
results, but that some care must be taken when using the latter results. 
The low Reynolds number results of Reference 1 are presented as a ratio 
of the low Reynolds number heating to the classical boundary layer heating 
and are correlated in terms of shock density ratio and Reynolds number. 

The results of the correlation are shown m Figure 7. The classical 
boundary layer solution used to make the low Reynolds number heat transfer 
results non--dimensiona 1 is based on modified Newtonian boundary layer edge 
properties. High density ratio (p /p^,) inviscid shock layer property 
results can differ considerably from modified Newtonian results. Thus 
the deviation of the low Reynolds number heat transfer from the so-called 
classical boundary layer result given m Figure 7 includes not only low 
Reynolds number effects, but also inviscid property effects. 

In addition to low Reynolds number effects, the effect of significant 
concentrations of CC^ in the ambient gas must be considered m the 
convective heat transfer calculations. Numerous theoretical and experi- 
mental studies of stagnation point convective heat transfer in N^-CO^ 
gas mixtures of interest in planetary entry have been made. The general 
result is that convective heat transfer increases as the percentage of 
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CO is increased. Predictions of the magnitude of the increase vary con- 

^ 23 

siderably, however. In this study the results of Scala and Gilbert have 

been used to account for the ambient gas composition. In their analysis, an 

extensive number of gas mixtures has been studied, and their results provide 

a reasonable upper bound on the effects of CO^ on convective heat transfer. 


The correlation results of Hoshizaki have been utilized to compute the 
stagnation point convective heat transfer, including the effecLs of dissocia- 
tion and ionization. The theory is applicable to either frozen or equilibrium 
chemistry boundary layers and the results are correlated in terms of shock 
properties and the inviscid flow stagnation point velocity gradient The re- 
sults of the analysis correlate closely to the existing experimental data. 

The blunt body wall has been assumed to be perfectly catalytic, and, for the 
flight velocities of interest m this study, the convective heat transfer is 
insensitive to the degree of non-equilibrium m the dissociated boundary layer. 
The correlation of the numerical solutions given in Reference 22 m terms of 
flight conditions is accurate to +1.57* m the flight velocity range of interest 
in this study The correlation formula in terms of stagnation point velocity 
gradient is given as 


q g (BTU/ft 2 sec) = 4.32 x 10 4 (- 


U 1.69 , T 0.4 % 

” . r dv . w . 

lO,OQO } L ds w U 'w K 500 ) 


(25) 


where U^, is the free-stream speed m ft/sec, dv/ds is the stagnation point 

velocity gradient in sec , p w is the density of air at the body surface m 
3 

slug/tt , p is the viscosity of air at the body surface in slug/sec-ft, and 


w 


T is the body surface temperature in K. A body surface temperature of 300 K 


is assumed for the present calculations. 


The heat transfer distribution around the body can be calculated by con- 
ventional boundary layer theory, as long as the appropriate boundary 

3 

layer edge conditions are used. For the flight conditions given in 
this study, the shock layer Reynolds number based on the blunt cone nose 
cap radius is approximately 9,000. At this Reynolds number, the boundary 
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layer thickness in the nose region is approximately 157* of the shock layer 
thickness.'*’ Since the boundary layer is a significant percentage of the 
shock layer thickness, the proper boundary layer edge conditions differ 
considerably from those given by inviscid flow body properties. The 
magnitude of the velocity at the edge of the boundary Layer is increased 
due to the shock curvature induced vorticity and can produce a significant 
increase m the convective heat transfer. The viscous and conduction 
effects due to velocity and temperature gradients in the shock layer, 
however, have a small effect on convective heat transfer at the wall and 
can be neglected. It should also be noted that the vehicle wall has been 
assumed to be cold, and the boundary layer displacement thickness and its 
interaction with the inviscid flow field is negligible. 

Based on a transition momentum thickness Reynolds number criterion of 250, 

the boundary layer will be laminar for the vehicles considered m this 

study. The laminar heat transfer distributions have been calculated by 

25 

Lees' local similarity solution. The Lees' method gives the ratio of 

the local heat transfer at a body station to the stagnation point value 

as a function of the local density-viscosity product and velocity at the 

26 3 

edge of the boundary layer. Cheng and Ferri have shown that, for 
spherical bodies at shock layer Reynolds numbers given m this study, the 
effect of vorticity on the laminar heat transfer distribution is approx- 
imately constant over the entire sphere cap. Thus, the inviscid body 
properties on the sphere-cap can be used in the method of Lees and the 
vorticity correction is obtained implicitly when the non-dimensionalized 
heat transfer distribution is multiplied by the low Reynolds number stag- 
nation point heat transfer result. 

Lees' analysis is also used to calculate the heat transfer distribution 
aft of the sphere-cone junction, but the inviscid properties at the edge 
of the boundary layer, rather than the vehicle wall , are used in the cal- 
culation. For the sphere cone cases, the edge of the boundary layer 
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properties in the conical region differ considerably from the mviscid 
body properties since the boundary layer edge extends a significant distance 
into the shock layer. 

Ideally, the boundary layer edge properties could be obtained by coupling 

3 

the boundary layer solution of Ferri et.al. with the detailed inviscid 
flow field results obtained m this study. Such a calculation requires 
considerable formulation and computer programming, however, and a 
simplified method that is estimated to deviate less than 5 % from the more 
exact analysis has been used. The boundary layer thickness in the conical 
region of the flow is assumed to be the same percentage of the shock 
layer thickness as the nose cap, namely 157=,. The inviscid properties at 
this station in the shock layer correspond almost identically to the 
inviscid flow field properties calculated in the cells adjacent to the 
body. The entropy region along the body induced by the blunt bow shock is 
completely contained between the body and the first cell and, hence, is 
swallowed by the boundary layer. Thus, the inviscid flow velocities at 
the boundary layer edges are more closely characterized by the conical 
shock entropy rather than the blunt nose normal shock entropy. 

2.3 RADIATIVE HEAT TRANSFER ANALYSIS 

From the results of the inviscid analysis of Section 2.1 the pressure and 
enthalpy of the gas throughout the inviscid portion of the flow field are 
obtained. From this information, three operations are carried out in 
determining the radiative heat transfer to the vehicle surface. The 
gas temperature and composition are first determined by the Aeronutronic 
free-energy-mimmization computer program at selected values of enthalpy 
and pressure that cover the range of these variables in the flow field. 

The resulting values of composition, temperature, and pressure are then 
used in the equilibrium radiation opacity program to obtain spectral 
absorption coefficients and radiation emission per unit volume. The 
radiation emission from each cell in the flow field is thus tabulated. 
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Finally the contribution of each cell to the radiation to selected points 
on the vehicle is computed and summed over all of the cells. In this step 
self -absorption is evaluated and appropriate corrections are applied. 
Radiation cooling corrections and radiation-viscous layer coupling correc- 
tions are then applied to obtain the final radiative heat transfer 
dis tribution . 


2.3.1 EQUILIBRIUM GAS COMPOSITION 

The free-energy minimization technique used at Aeronutronic to determine 

equilibrium gas properties has been described by Oliver, Stephanou, and 
27 

Baier, and is now a commonly used method. The detailed thermodynamic 
properties of the major species used have been tabulated by Hildenbrand, 
and have been maintained current with the JANAF Tables. For planetary 


28 


atmosphere calculations, the program includes the species C, N, 0, A, C^, 

N 2 , 0 2 , C0 2 , CN, CO, NO, C 3 , e", 0", 0 2 ", A*", C + , 0 + , N + , C 2 + , N 2 + , 0 2 + , 

CN + , C0 + , N0 + , and C (solid). Multiply ionized atoms are also allowed, 

but are not found in significant concentration in the present calculations. 

The effects of vibrational anharmomcity , centrifugal stretching, and 

rotation-vibration coupling are included in the partition functions of 

the molecules. The known low-lying excited states of all species are also 

27 

included. The program uses the free-energy minimization principle to 
determine the temperature and composition for any prescribed gas enthalpy 
and pressure. 


The heat of formation of CN is listed m the JANAF Tables as 109 Real/ 
mole, corresponding to a dissociation energy of 7.5 ev. In the present 
calculations a value of 91 Kcal/mole was used, which corresponds to a 
dissociation energy of 8.2 ev. The latter value is consistent with the 
shock tube measurements conducted at Ames, and is consistent with the 

shock tube measurements of CN red emission carried out by Thomas and 

. 29 29 

Menard at JPL. Although some CN violet measurements suggest a lower 
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value for the dissociation energy, the higher value has been selected 
for the present calculations. 

The results of the Aeronutronic free-energy-minimization program have 
been compared with those of JPL and of Ames and are found to be in 
satisfactory agreement. Species concentrations differ by as much as 10 
to 207o m some cases, due to differences in the thermal properties of the 
species that cannot be resolved on the basis of presently available 
information . 


2.3.2 EQUILIBRIUM RADIATION EMISSION 

The equilibrium radiation emission program used for the present study 

30 31 

is described m detail in Aeronutronic reports by Bauer and Main. 5 

The frequency-dependent absorption coefficient of CO„-N mixtures is 

- 1 1 -1 

computed over the spectral range 1000 cm to 67,000 cm for 
equilibrium temperatures from 1000°K to 10,000°K. The volume radiation 
emission is calculated from the absorption coefficient and the falackbody 
radiation function. 


An extensive search for all known data on equilibrium radiation from 
N^-CO^ mixtures has been carried out specifically for the Martian entry 
problem. The principal sources of radiation are the electronic band 
systems of the diatomic and polyatomic molecules, the vibration-rotation 
bands, and the free-free and free-bound continuum. The techniques and 
data used m the program are summarized below. 

a. Molecular Electronic Transitions . This mechanism contributes 
most of the radiation emission for the temperatures and pressures 
associated with significant radiative heat transfer during Martian 
entry. Here a diatomic or polyatomic molecule in an electronic state 
E 1 , vibrational state v', and rotational state J ( , emits a quantum of 
radiation to go to a state (E", v", J") . The presence of large numbers 
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of different v- and J-values gives rise to the characteristic structure 
of these transitions or "band systems" which are usuaLly quasi-cont mua 
at the temperatures considered here. 


The model used for the band systems of diatomic molecules is called the 
just overlapping rotational line ("JORL") band model or the smeared 


rotational line band model. 


32 


This band model assumes the Born- 


Oppenheimer separation of electronic and nuclear motions in the molecuLe 

and the perturbation treatment of the coupling of its vibrational and 

33 

rotational degrees of freedom used m spectroscopic work. A profile 
is assigned to each band of this system corresponding to that of the Q- 
branch of the band and containing all of the integrated strength of the 
band. While this model does not give the correct profile of each band, 
it does give an accurate representation of the total absorption and 
emission of the band system. A simple extension of this procedure is 
-applied to polyatomic molecules, e.g., CC^ • 


33 


To compute the equilibrium radiation emission from each band system, 
the Franck-Condon factors and the oscillator strength are required in 
addition to the molecular spectroscopic constants. The Franck-Condon 

34 

factors are computed using the method of Nichols for those band systems 

35 

which have not been computed using more rigorous methods. The electronic 

oscillator strengths are based on published measurements or calculations 

wherever possible. In the absence of previously published f numbers, 

values have been estimated using published values for similar allowed 

36 

transitions in isoelectronic molecules. In other cases the qualitative 

37 

rules of Muiliken and Rieke are of assistance in arriving at estimates. 
References and values for the Franck-Condon factor arrays and the system 
oscillator strengths used in the radiation program are Listed in Table 1. 

b. Vibration-Rotation Bands . The vibration rotation bands of NO, 

CO, and CN contribute significantly to the radiation emission at tempera- 
tures below 3,000°K. These bands are included m the radiation program 
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TABLE 1 

DATA FOR ELECTRONIC BAND 
MOLECULES IN C-N- 


1. 


2 . 


3. 


4. 


6 . 


7. 


Band System 

(first positive) 

N 2 (C 3 n u -B 3 n g ) 
(second positive) 

N 2 (B ' 3 V B \ ) 

("Y" system) 

+, 2„+ 2--K 

L(B Z -X 2 ) 

2 u g 

(first negative) 

^(A^-xV) 

(Meinel) 

N+( C 2 2+-xV) 

(second negative) 

< (d V* 2 v 


SYSTEMS OF DIATOMIC 
•0 MIXTURES 

System Absorption f- 
number and References 


0.0028 


57 j 5S 


0.043 59 ’ 60 

0.002 (est) 

36 ,60,61 
0.035 1 5 

0.007 (est) 

0.10 (est) 

0.005 (est) 



(Janin-D* Incan) 


8. 

(Schumann-Runge) 

0.17-(0,34xl0 4 )/^ 

_ L 1 

S 200000 cm ) 

9. 

N0(B 2 n r -X 2 TT r ) 

(beta) 

59 

0.0053 

10. 

2-4-2 

NO(A 2 -X n r ) 
(gamma) 

0.OO25 59,62 

11. 

NO(C 2 n-X 2 n r ) 
(del ta) 

0.017 62 

12. 

R0(D 2 2 + -X 2 TT r ) 

(epsilon) 

O.OH 62 ' 63 
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TABLE 1 (continued) 



Band Systems 

System Absorption £- 
number and References 

13. 

N0(E 2 £ + -X 2 rr r ) 

0.01 (est) 


(gamma prime) 


14. 

NO(B ,2 A -X 2 tt ) 
(beta prime) 

0.01 (est) 

15. 

NO(G 2 Z"-X 2 Tt ) 

0.01 (est) 


(Lagerquist-Miescher) 

16. 

N0(C 2 tt-A 2 E + ) 

0.05 (est) 


(Heath) 


17. 

no(d 2 £ + -a 2 s + ) 

0.05 (est) 


(Feast 1) 


18. 

K0(E 2 S + -A 2 tt) 

0.05 (est) 


(Feast 2) 


19. 

N0(E 2 I + -C 2 n) 

0.05 (est) 

20. 

no(e 2 e + -d 2 " + ) 

0.05 (est) 


(Feast-Heath) 


21. 

2 2 

N0(B 1 a -B tt ) 
l r 

0.001 (est) 


(Ogawa 1) 


22. 

NO(G 2 £”-B 2 TT r ) 

0.001 (est) 

23. 

NO (b^£ -a^n^) 

0.001 (est) 


(Ogawa 2) 

' 

24. 

N0 + (A 1 tr-X'*‘2' f ) 

0.05 (est) 


(Miescher-Baer) 


25. 

CN(B 2 S ! '-X 2 E + ) 

0.027 M > 65 


("violet" system) 
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TABLE 1 (continued) 



Band Systems 

System Absorption f- 
number and References 

26. 

CN(A 2 Tr^-X 2 E + ) 

0.007 66 


("red" system) 


27. 

CN(D 2 Tr^-X 2 E^) 

0.005 (est) 

28. 

CN(D 2 TT.-A 2 n.) 
l i 

0.005 (est) 

29. 

CN(F 2 A -A 2 tt.) 
' r l 

0.01 (est) 

30. 

CN(E E -A rr ) 

0.01 (est) 

31. 

CN(J 2 A.-A 2 tt. ) 
ii 

0.01 (est) 

32. 

CK(E 2 E + -x 2 E + ) 

0.10 (est) 

33. 

CN(H 2 n r -X 2 E + ) 

0.05 (est) 

34. 

C 2 (A 3 w g -X' 3 n u ) 

0.034 67 ' 68 


(Swan) 


35. 

C 2 (A' 3 2--X' 3 n u ) 

(Ballik-Ramsay) 

0.005 (est) 

36. 

37. 

C 2 (b 1 V xlE g ) 

(Phillips) 

69 

0.0038 

37. 

c 2 (B 3 ng-X' 3, i u ) 

(Fox-Herzberg) 

0.13 (est) 

38. 

_ . ,l v + 1 V + N 
C 2 (d E u' x V 

0.10 (est) 


(Mul liken) 


39. 

c 2 (e E g -b tt u ) 
(Freymark) 

0.10 (est) 

40. 

C 2 (c*tt -b^rr^) 0.005 (est) 

(Deslandres-D'Azambuja) 
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TABLE 1 (continued) 



System Absorption f- 


Band Systems 

number and References 

41. 

CO(A 1 tt-X^T" } ') 
(fourth positive) 

0.15 63 

42. 

CO + (A 2 tt.~X 2 £ + ) 

0.0022 70 


( "comet- tail" system) 

43. 

co + (B 2 r + -x 2 s + ) 

0.01 (est) 


(first negative) 

. 

44. 

i _^+, t , 2 v + .2 . 

CO (B S -A n i ) 

(Baldet- Johns on) 

0.001 (est) 

45. 

C0(B L 2 + -A 1 n) 

(Angstrom) 

0.005 (est) 

46. 

C0(C I - + -A I tt) 

(herzberg) 

0.005 (est) 

47. 

CO(a ,3 Z + -a 3 Ti) 

0.015 (est) 

- 

(Asundi) 


48. 

3 + 3 

C0(b £ -a n) 
(Third positive) 

71 

0.009 

49. 

C0(d 3 A -a 3 n) 
(Triplet) 

0.015 (est) 

50. 

CO(e 3 l”-a 3 n) 

(Herman) 

0.015 (est) 

51. 

4-9 9 

0:(A‘ tt -X tt ) 

2 u g 

(Second negative) 

0.005 (est) 

s 52. 

°+( b 4 r-a\> 

(First Negative) 

0.001 (est) 
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by the following simple treatment. It has been shown theoretically that 

the integrated absorption coefficient of the fundamental system of a diatomic 

molecule is independent of temperature when induced emission is taken into 

account. This has been demonstrated to be true m experimental measurements 

38 39 

for the NO and CO systems. 5 Using these data, and an estimate for CN, 
the average absorption coefficients over the spectral range 1,000 to 3,000 
cm 1 , p, (cm 1 ), are 13.6 y nq P/T, 39.4Y CC| P/T, and 24.6 y cn P/T, where Y^ is the 
mole fraction of the species, i, P is the total gas pressure (atm) and T is 
the gas temperature (°K) . 

The integrated intensity of the first overtone system of diatomic molecules 

38 

has been shown to increase with increasing temperature theoretically and 

39 40 

experimentally for NO and CO . * We take the average absorption coefficients 
of NO, CO and CN in the spectral range 3,000 to 5,000 cm ^ to be 

i , 

p (NO) = (0.000228T + 0.182) Y nq P/T 
p'(CO) = (0.000195T + 0.0778) Y P/T 

p'(CN) = 0 .346[ 1- exp (-20 69 hc/kT)]" 2 [l-exp(-4138 hc/kTJY CN P/T 

The value for CN is an average of the measured values for NO and CO at 
2,000°K, and the theoretical expression of Reference 38 is applied. 

The absorption coefficients defined above are satisfactory only in computing 
optically thin emission Each vibration- rotation band actually consists of 
50 to 100 rotational lines. All of the lines cover a spectral range that is 
5 to 10% of the nominal wavelength of the band. However each rotational line 
covers a very narrow portion of the spectrum at the conditions of temperature 
and density associated with these flight calculations. The Doppler broadened 

i, g 

linewidth is (2nkT/M) 2 A and is approximately 3x10 cps at the stagnation 
point for the present flight case. Collision broadening yields a linewidth 

g 

of about 10 cps, contributing less than Doppler broadening. Since the 
Doppler linewidth represents only 10 ^ of the nominal frequency of the 
vibration-rotation bands, these lines are so sharp that self -absorption is 
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significant even when the absorption coefficient defined above is much less 

-3 

(about 10 ) than the reciprocal of the thickness of the shock layer. It 

can be seen by inspecting Tables 6 and 7 that the value of p ' at 2000 and 
4000 cm ^ approaches, but does not exceed the reciprocal of the shock layer 
thickness. Thus the infrared emission from the NO, CO, and CN bands 
approaches being self-absorbed. However, since the magnitude of this emission 
contributes less than 10 % of the total radiation to the body, any errors 
introduced in the total radiation by the omission of self-absorption correc- 
tions in the* vibration-rotation bands must be considerably less than 107 o . 


Mo provision was made to include the vibration-rotation bands of CO^* For 
the equilibrium case, the concentration of CO^ is so small that the emission 
is negligible. For the non-equilibrium cases, the CO^ volume emission is 
appreciable immediately behind the shock prior to CC >2 dissociation, but when 
integrated over the shock layer it is found to be negligible in comparison 
with the integrated emission due to the CN red and CN violet bands. 

c. Continuum Emission . Continuum emission due to free-free transi- 
tions is computed using the classical Kramers formula for the linear 
spectral absorptioji coefficient, 

2 

Y y P 

= 0 . 738xl0 21 z^ f f | 5/2 tl- e xp(-hcw/kT)J 

T 


which includes allowance for stimulated emission. The species is taken 

here to be either the totality of all positively charged singly ionized 

2 

atoms and molecules for which we take z = 1 , or the totality of all 

eff 2 

electrically neutral atoms and molecules, for which we take z“ f = 0.01 

42 

following Keck, Allen and Taylor. The continuum emission from free- 

bound transitions is determined by using the hydrogenic cross sections 

43 

for recombination as computed by Bates. The continuum emission is 
appreciable for mixtures only when the equilibrium temperature 

exceeds 8,000 to 10,000°K. 
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2.3.3 SELF- ABSORPTION OF EQUILIBRIUM RADIATION 

For equilibrium calculations, self-absorption corrections to the radiative 
heat transfer are made in two different ways. For bands such as 00(4+) 
that are very optically thick a blackbody upper limit is set at each wave- 
length by utilizing the temperature of the shock layer gas at a normal 
distance of one absorption length from the body surface. The absorption 
length for each wavelength is printed out directly from the equilibrium 
radiation program. The blackbody limits at each wavelength are then 
summed over all wavelengths. 


For intermediate optical depth the following model is used to make self- 
absorption corrections to the radiative heat transfer. The spectral 
emission (neglecting self- absorption) is computed and integrated progressively 
from the body surface in each direction defined m Paragraph 2.3.6. 
Simultaneously the non-equilibrium absorption optical depth is integrated. 

At a point where the absorption optical depth reaches a value of one, the 
integration of the optically thin emission is terminated; the total emission 
thus accumulated closely approximates the net spectral radiative flux with 
self-absorption . 


2.3.4 RADIATION- INVISCID FLOW COUPLING 

This is negligible in the equilibrium sphere-cone case and is therefore 
not included. The radiation cooling of the inviscid fluid by the time 
it reaches the boundary layer edge adjacent to the stagnation point 
reduces the enthalpy by about 1 % and reduces the total radiation emission 
by less than 1%. 
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2.3.5 RADIATION-VISCOUS LAYER COUPLING 

Radiation-viscous layer coupling is included by subtracting from the total 
radiation a contribution from a uniform slab of thickness equal to the 
thermal boundary layer thickness, defined by a 57, reduction in enthalpy 
below the value just outside the boundary layer. The radiative properties 
of the slab were taken to be those in the cell nearest the body. The addi- 
tional effect of radiation absorption by the boundary layer gas is neglible, 
because the enthalpy flux through the boundary layer gas is much greater 
than the radiation absorbed in the boundary layer. 

2.3 6 GEOMETRIC INTEGRATION 

The radiative properties of the shock layer gas determined from the 
equilibrium radiation emission program are used in a geometrical integra- 
tion program to compute the radiative flux to various points on the 
vehicle surface. For a given point on the body, the shock layer is 
divided into slabs as shown in Figure 8. These slabs have their normal 
coincident with the normal to the body through the body point for which 
the radiative heat transfer is being computed. In the present computation 
three slabs of equal thickness were used between the body and the shock. 

Each slab is divided into rings that subtend equal solid angles from 
the body point. Each ring in each slab thus contributes equally to the 
flux to the body point if the radiative properties throughout the slabs 
are uniform. However, because of the shock curvature, the outermost 
slab is only partially filled with r adiating gas. Furthermore, the shock 
layer properties are not uniform These two effects are included by 
summing the radiative contributions from segments of each ring. The 
contribution from each segment is based on an average radiating volume 
determined from the geometry of the slabs and the shock layer, and an 
average radiation emission per unit volume from the segment. The 
summation of these contributions gives the total radiative flux to the 
surface. In the present computations each slab was divided into five 
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FIGURE 8. 


GEOMETRY USED IK CARRYING OUT THE VOLUME 
INTEGRATION OF SHOCK LAYER RADIATIVE HEAT 
TRANSFER TO SELECTED BODY LOCATIONS 


- 39 - 



rings, the first two nearest the normal subtending one half the solid 
angle subtended by the other three. The angles used to compute average 
radiative properties were 14.5°, 35.25°, 50.75°, 67.75°, and 82.75°. 

2.3.7 COLLISION LIMITING 

It has been found from the non-equilibrium studies of Section 3.2 that 
the equilibrium CN violet emission may be reduced by as much as a factor 
of 2 for the conditions of the present calculations. The other significant 
sources of radiation do not appear to be affected. However because of 
the uncertainty m this factor, and because it has not been corroborated 
by experimental measurements, no correction for collision limiting is made 
to the equilibrium radiation calculations. 
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SECTION 3 


NON- EQUILIBRIUM FLOW 


The flow field analysis used for non-equilibriUm conditions is much the 
same as that used for equilibrium conditions, described m Section 2. 

The inviscid properties are computed assuming a thin viscous layer, and 
are then corrected for coupling between the inviscid and viscous regions. 
The convective heating and radiative heating are then calculated with 
the use of, the uncoupled inviscid results, corrected for radiation 
cooling and viscous-inviscid coupling. 

3.1 INVISCID ANALYSIS 

In the inviscid region the flow equations are decoupled from the reaction 
rate equations by including an ideal gas equation of state in the flow 
equations. This separation permits an effective Y to be used which is 
determined through the results of an integration of the reaction rate 
equations along the stagnation streamline. The variation of Y along the 
stagnation streamline is found to be representative of that along other 
normals between the body and the shock throughout the subsonic flow 
region in the shock layer. In this way the equation of state differs 
little from that obtained from a rigorous solution of the flow equations 
coupled with the reaction rate equations. In the present calculations 
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the y variation noted above was approximated by a constant value since 
the variation in y is not Large across the shock layer. 

The results of integrating the flow equations with the approximate 
equation of state are used only to provide the pressure field and the 
streamline positions. The remaining gas properties are determined by 
solving the reaction rate equations along with Euler's equation, and the 
equations of state, on selected streamlines. The errors introduced by 
this procedure are small since the pressure field and streamline posi- 
tions are insensitive to non- equilibrium chemistry effects. 

In the constant y flow field solution, the conservation Equations (1)- 
(5) are solved by the same method used for equilibrium flow (Section 2.1). 
The thermal and caloric equations of state are different, however, and 
for constant y combine to give 

P = (Y-l)pe (26) 

which is used in place of Equation (8). 

The streamline positions within the shock layer are determined by locating 
points of constant stream function, \jj, which by definition are points on 
a given streamline. Lines of constant stream function were calculated 
from an integral solution of the continuity equation utilizing the density 
and velocity profiles obtained from the mviscid shock layer solution. 

The non-equilibrium chemistry calculation is initiated at the bow shock 
wave. A.cross the bow shock the rotation and the vibration internal energy 
modes are assumed to be fully equilibrated with translation, while the 
chemical species are assumed to be frozen at the free stream conditions. 

The y utilised in the inviscid flow field solution corresponds to a 707, 

$ 2 > 307c CO 2 mixture with equilibrium internal energy modes and frozen free- 
stream chemistry. Thus, the initial conditions for p, p, T, q, Y^ and E^ 
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behind the bow shock to start the non-equilibrium species calculation are 

given by the inviscid solution. With the pressure and streamline field 

given, the remaining five unknown shock layer properties, p, T, q, V j , and 

E are determined from the simultaneous solution of the following five 
v . 

equations: 


Euler Equation 


dq _ -k dg 
dz pq dz 


( 27 ) 


Integral Energy Equation and the Caloric Equation of State 
s £ 


T = 


k <±12 : i£i ( W >h ° - chVyns 


? V 1 - 5+ V 


Thermal Equation of State 


P * 


(T £ y ) 


( 28 ) 


(29) 


Species Conservation Equation 

dy . s Q x . 
q dz x ~\ p 

where 


( 30 ) 


Q v . 1 v i-r 

-ii . (i <c P ) [l 


j3 . . = v 1 - v . 

ij xar io? 


, P s B 
z-(cp) 1 11 Y 3 
K. o-=l 

i 
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Equilibrium Vibration Energy Equation 


~ ~ 1 u 

E = RT[u(e U -l)" A + — (31) 

v — 

u 

0 

where u = — and the constants u and 0 are given in Table 2. 

T o v ° 

3.2 NON- EQUILIBRIUM PROPERTIES OF PLANETARY ATMOSPHERES 
An extensive compilation of non- equilibrium processes in CO^, A 
mixtures is included in the Aeronutronic non- equilibrium computer program, 
specifically for Martian entry calculations. The program includes the 
processes normally included in all such programs: molecular rotation and 

vibration, dissociation, ionization, and atom, ion, and charge transfer. 

In addition, careful consideration has been given to the most important 
process for computing radiative heat transfer, that of non-equilibrium 
excited state production. 

The non-equilibrium properties of the flow field are obtained from a 
streamtube integration program which utilizes the pressure field and the 
streamline locations determined from the flow field integration (Section 
3.1). The streamtube program integrates the energy and Euler equations, 
and the species and excited state rate equations to obtain the chemical 
composition, the gas density, and the temperature. 

The details of the program and the non-equilibrium data used in it are 
summarized below. The species included in the program are C, N, 0, , 

N 2 ’ ° 2 * C 0 2 ’ C0 > CN > K0 > C+s N+ > 0+5 C 2 + > N 2 +j ° 2 + ’ 00+5 CN+j N0+ ’ e ' 5 
CN*, N 2 *, and CO*. The number of reactions included is 31 for pure 

93 for N 2 -CO 2 mixtures, and 27 for air. 



TABLE 2 

' MOLECULAR VIBRATION AND DISSOCIATION CONSTANTS 


Molecule 

o 

r\ ^ XT 

u 

D , 

kcal/ mole 

Refs . 

9^ K 

— *CT“ 

J 


CO 2 = 1950 

0^ - 966 

0 3 = 3445 0 

126 

(to CO + 0) 

72,73 

CO 

3067 

.029 


256 

33 

N 

3337 

.032 


225 

33 

2 

0 

2223 

.039 


118 

33 

2 

n 

2322 

.032 


140 

33,74 

2 

CN 

2924 

.030 


172 

33,75 

NO 

2688 

.037 


150 

33 
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The reaction rate constants are based on an extensive survey of available 
experimental data. However there are a number of reactions which con- 
tribute significantly to the non-equilibrium properties of the flow, for 
which there is no data available. For these reactions the best methods 
available have been used to estimate rate constants. The basis for these 
estimates and the values used in the program are discussed below. The 
reactions are divided into nine groups, which are identified as rotation 
relaxation, vibration relaxation, molecular dissociation and atom recombina- 
tion, atom exchange, neutral atom transfer, electronic excitation, ion 
production and loss, and ion and charge transfer. The data used m the 
non-equilibrium program are given here m detail because they are not yet 
available elsewhere. 

3.2.1 REACTION RATE DATA 

a. Rotation Relaxation . Rotation is assumed to be fully equilibrated 
with translation, as has been observed m low density shock thickness 

44 45 

studies in CO^ and 0^ • The rotational contribution to the heat 
capacity is RT for all molecules, including CO^ . 

b. Vibration Relaxation. For normal shocks above V_ = 4.5 km/sec, 

CO and vibrational modes are fully relaxed within the translation- 
^ ^ 44 45 

rotation shock. ’ Vibration relaxation times for and CO are within 
a factor of 5 of that for 0 2 above 8000°K. Furthermore vibration is 
in resonance with the mode of CO^ and will relax even faster in CO^ - ^ 
mixtures. Thus for normal shock speeds above 4.5 km/sec vibration is 
fully relaxed within or very close to the translation-rotation shock, and 
can be assumed to be fully equilibrated with translation. This applies 
as well to oblique shocks with a normal component of speed greater than 
4,5 km/ sec. 

The equilibrium vibrational energy is obtained from the spectroscopically 
determined vibrational constants by = RT[u(e U -1) ^ + u q /u] where u = 
G v /T, (see Reference 46). Values for 9^ and are given in Table 2. 
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Only the vibrations of the neutral species C0 2 , CO, N 2 , 0 2 , C^, CN, and 
NO are considered to contribute to the gas heat capacity. The special 
case of CO^j having 4 vibrational modes (2 modes at 1 frequency, v 2 )» is 
treated by uncoupled superposition of each mode, neglecting the effects 
of anharmonicity (u^) since the constants \ and y are quite small for CO^ . 
Because the \> 2 (bending) mode is doubly degenerate it contributes = 

-1 ' >2 
2RTtu(e U - 1) ] to the heat capacity. 

c. Molecular Dissociation and Atom Recombination . Dissociation of 
the neutral molecules is included, since these cause the major change 

in gas temperature. In addition dissociation of the excited electronic 
states, CN*, CO*, and N^* is included, since dissociation strongly affects 
the concentration of these species. Certain rates are estimated because 
of the lack of experimental data. This is done by assuming that molecules 
with paired electrons (0 9 , N 2 , C 2 , CO, and C0 2 ) are less efficient in 
causing vibrational energy transfer near the dissociation limit than are 
molecules with unpaired electrons (NO and CN) and free atoms (C, N, and 
0) . The former group is denoted by the latter by M 2> These rate 
constants are defined for a reaction of the form A + A + M -* 

A 2 + M by d[A 2 l/dt = - 1/2 d[A]/dt = k f l_A] [MJ . The reactions and the 
rate constants used m the Aeronutronic computer program are given as 1 
through 14 in Table 3. 

d. Atom Exchange . The atom exchange reactions are generally slow 
and exhibit a high activation energy even for exothermic reaction rates. 
Because of this, estimates cannot easily be made. Thus only those reac- 
tions for which data exist are included. These rate constants are de- 
fined for a reaction of the form AB + AB -* A 9 + B 2 by d[A 9 J/dt = 

- 1/2 dfAB]/dt = k c [AB]^. The reaction and rate constant used is given 
f 

as 15 in Table 3. 

e. Neutral Atom Transfer . The atom transfer reactions are generally 
fast and exhibit small (if any) activation energy for the exothermic reac- 
tion rates. Only a feu of these reaction rates have been measured. All 
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TABLE 3 

REACTION RATE SYSTEM FOR C-N-0 MIXTURES 



Reaction 

Forward Rate Constant, 
cm /mole-sec 

Ref 

Equilibrium Constant 

I. 

C0 2 + CO -h ° + M l 

0.116xl0 12 T°‘ 5 exp(-34,300/T) 

76 

0 . 41xlO S T _ 1 * 45 exp (- 63 , 000/T) 

2 . 

CO 2 + M 2 - CO + 0 + M 2 

3 k l 

est 

K 1 

3. 

0 + 0 + M l - 0 2 + M 

0.80x10 20 T~ 1 ' 5 

77 

0.64x10 5 T lr °exp(62,419/T) 

4. 

o + o + m 2 - o, + m 2 

2.75 k 3 

77 

K 3 

5. 

N + N + - N 2 + M 

0.22x10 17 t"°' 5 

78 

0 . 91x10 3 t" 1 ‘ °exp( 107, 365/T) 

6. 

N + N + M 2 - N 2 + M 2 

22 -1 5 
0.22x10 T 

78 

K 5 

7. 

N + 0 + M 1 - NO +• 

21 -1.5 
0.40x10 T 

est 

0.39xl0 2 T' 0 ' 5 exp(71,673/T) 

8. 

N + 0 + - NO + M 

7.5 k 

79 

K ? 

9. 

C + 0 + M - CO 4- M 

O.SOxlO 19 !' 1 * 0 

80 

0.23xl0" 3 T°‘ 5 exp(131, 119/T) 

10. 

C + o + m 2 - CO + m 2 

10 k 9 

est 

K 9 

11. 

C + N +■ M x - CN -h M 1 

0.25xl0 15 

est 

0.27exp(92,400/T) 

12. 

C + 11 F - CN f M 2 

6 k ll 

est 

K 11 

13. 

C + C + M x - C, + M x 

0.20xl0 17 T"°' 5 

est 

0 .23exp(66, 691/T) 

14. 

c + c + m 2 - c 2 + m 2 

3 k 13 

est 

*13 

15. 

NO +■ NO - N 2 + 0 9 

0.23xl0 24 T" 2 * 5 exp(-43,000/T) 

77 

K 3 K 5 /K 7 2 

16. 

co +■ o 2 - co 2 + 0 

0.10xl0 12 T°‘ 5 

est 

1/K 1 K 3 



TABLE 3 (continued) 


Reaction 
17. C 2 + 0 -* CO 



Forward Rate Constant, 
cm /mole-sec 

Ref 

Equilibrium Constant 

+ 

c 

0.50xl0 l2 T°* 5 

est 

V K 13 

18. 0 + C -* CO 

+ 

0 

0.50x10 12 T°' 5 

est 

V K 3 

19 . CN + 0 -* CO 

+ N 

0.10x10 13 T °' 5 

est 

K 9 /k 11 

20. NO + C - CO 

+ N 

0.25x10 12 T°' 5 

est 

V K 7 

21. 0 2 + N - NO 

+ 

0 

0.13xl0 11 T 1 ‘°exp(-3560/T) 

77 

K 7 /K 3 

22. CN + N - N £ 

+ 

c 

0.125x10 12 T°‘ 5 

est 

K 5 /K U 

23. NO +• N N 2 

+ 

0 

14 

0.16x10 

77 

K 5 /K 7 

24. NO 4- C - CN 

+ 

0 

0.125xl0 12 T°' 5 

est 

K ll /K 7 

25. C 2 + N - CN 

+ 

c 

0.50x10 12 T 0,5 

est 

K ll /K 13 

26. N + 0 - NO + 

+ 

e 

0.75xl0 14 T"°‘ 5 exp(-31,900/T) 

24 

0 . 66x10" V ' °exp(-35 , 993/T) 

27. 0 + 0 - 0 ^ 

+ 

e 

0.50xl0 l4 T'°' 5 exp(-82,300/T) 

24 

0.38x10' 13 T 2 ’ 5 exp(- 76, 387/T) 

28. N +■ N - N 2 

+ 

e 

0.15xl0 15 T~°' 5 exp(-67,500/T) 

24 

0.3lxl0" 6 T 1,O exp(-74,O51/T) 

+■ 

29 . C + N - CN 

+ 

e 

O.8Oxl0 l4 T"°‘ 5 exp(-69,6OO/T) 

est 

0.42xl0“ 6 T 1 ’°exp(-84,122/T) 

30. C + 0 - CO + 

+ 

e 

0.27xl0 14 T"°- 5 exp(-33,600/T) 

es t 

0.32xlO' 13 T 2 ' 5 exp(-29,370/T) 

31. C + C - C 2 + 

+ 

e 

0. 80xl0 li+ T _ °‘ 5 exp(- 68 , 600/T) 

est 

0.30xl0" 6 T 1,0 exp(-75,950/T) 

32. 0 + + NO - 0, 

+ 

7 

+ N 

0.25x10 12 T°* 5 

est 

0.29xl0" 6 T 1,5 exp(9437/T) 



TABLE 3 (continued) 



Forward Rate Constant, 
f’m / mole- sec 

Ref 

Eauilibrium Constant 

Reaction . — 

33. 0 + + N„ NO + + N 

0 . 50 x 10 12 T °' 5 

est 

K 26 K 32 /K 23 K 27 

2 

4 " 4 " 

34. N + N 2 N 2 + N 

0 . 25 x 10 12 T °* 5 

est 

0 . 53xlO^T~^ '° exp (11, 113 /T) 

35. 0 + + 0 2 - 0 + + 0 

0 . 50 x 10 12 T °* 5 

est 

K 21 K 32 

36. N + + NO - N 2 + + 0 

0 . 25 x 10 12 T °‘ 5 

est 

K 23 /K 34 

37. 0 + + NO - NO + + 0 

0 . 25 x 10 12 T °* 5 

est 

K 23 K 33 

38. N + + 0 2 - 0 2 + + N 

0 . 50 x 10 12 T °' 5 

est 

K 5 K 27' /K 3 K 28 K 34 

39 . C + 0. + -* CO + + 0 

0 . 25 x 10 L 2 T °* 5 

est 

K 30 /K 27 

40. C + + 0 - CO + + 0 

0 . 50 x 10 12 T °' 5 

est 

0 . 45 x 10 + 5 T _ 1 ‘° exp ( 24 , 943 / T ) 

4 

4 * "h 

41. 0 + CO - CO + o 

0 . 25 x 10 12 T °' 5 

est 

K 18 /K 35 K 39 

4- 4* 

42. 0 + CN -» CO + N 

0 . 25 x 10 I 2 T °’ 5 

est 

K 19 /K 41 

43. 0 + CN + -» CO + + N 

0 . 125 x 10 12 T °‘ 5 

est 

K 30 /K 29 

44. C + + NO - CO + + N 

0 . 25 x 10 12 T °‘ 5 

est 

k 40 / k 2 i 

45 . N + CO + - N0 + + C 

0 . 125 x 10 12 T °* 5 

est 

K 26 /K 30 

46. 0 + + C 9 -» CO + + C 

4- “1" 

47. C + CO -> CO + C 

0 . 50 x 10 12 T ° 

est 

K 17 /K 41 

0 . 125 x 10 12 T ° 

est 

K 20^ K 44 

48. 0 + c + “* CO + + C 

0 . 25 x 10 12 T ° ‘ 5 

est 

K 3l^ K 32 



TABLE 3 (continued) 


Forwaid Rate Constant, 


Reaction 

cm /mole-sec 

Ref 

Equilibrium Constant 

49 . 0 + C 2 + - CO + C* - 

0 . 25 x 10 12 T -°’ 5 

est 

o 

CO 

CO 

r— < 

50. C + CN + - C + + N 

0 . 125 x 10 12 T ° 

est 

K 31 /K 29 

51. N + C 2 + - CN + C + 

0 . 25 x 10 12 T °' 5 

est 

K 49^ K 19 

52. C + C 2 + - C 2 + C + 

0 . 25 x 10 12 T °* 5 

est 

K 51 /K 25 

53.0+ CN + - CO + N + 

0 . 125 x 10 12 T °' 5 

est 

K 18 K 27' /K 29 K 38 

54. 0 + CN + - N0 + + C 

0 . 125 x 10 12 T °‘ 5 

est 

K 26 /K 29 

55 . 0 + CN + -» NO + C + 

0 . 125 x 10 12 T °' 5 

est 

K 21 K 43^ K 40 

56. N + + CN -* CN 4 " + N 

0 . 25 x 10 I 2 T °‘ 5 

est 

K 19 /K 53 

4- 4- + C 

57. W + CN - N„ 

0 . 125 x 10 12 T °‘ 5 

est 

K 28^ k 2 9 

4- Z + 

58. N + CN - N 2 + C 

0 . 125 x 10 12 T °‘ 5 

est 

K 5 K 50 K 52' /K 13 

59. N + + C 2 - CN f + C 

0 . 50 x 10 12 T °' 5 

est 

K 25 K 56 

60. C + CN + - CN + C + 

0 . 125 x 10 12 T °* 5 

est 

k 58 / k 22 

61. N + 0^ “* NO 4 ' + 0 

0 . 25 x 10 12 T °’ 5 

est 

K 26 /K 27 

62. N + + 0, -* N0 + + 0 

0 . 50 x 10 12 T °' 5 

est 

K 38 K 61 

63. N + + NO - NO 4 " + N 

0 . 25 x 10 12 T °‘ 5 

est 

K 62 /K 21 

64. 0 + N * - N0 + + N 

0 . 25 x 10 12 T °* 5 

est 

K 26 /K 28 

65. 0 + + CN NO 4 " + C 

0 . 25 x 10 12 T 0,5 

est 

K 37 /K 24 



TABLE 3 (continued) 


1 

Ln 

N> 

I 


Reaction 

66 . + CO - NO + C 

67. N* 4* CN - + C 

68 . C + 0 2 + - CO + 0 + 

69. N + "* N 2 * + N 

70. N + CN - N 2 * + C 

71. N + NO -* N 2 * +■ 0 

72. C + CO - CO* + C 

73 . 0 + CO -* CO* + 0 

74. C + 0 2 - CO* + 0 

75. 0 + C 9 - CO* + C 

76. C + CN - CO* + N 

77 . 0 + CN - CO* + N 

78. C + C0 2 -* CO* +• CO 

79 . N 2 * + CO - CO* + N 2 

80 . C + 0 + M x c0 * + M l 

81. C + 0 + M 2 - CO* -6 M 2 


Forward Rate Constant, 
cm / mole-sec 

Ref 

Equilibrium Constant 

0.25x10 12 T°’ 5 

est 

K 62 /K L8 

0.25x10 12 T°- 5 

est 

K 22 /K 34 

0.25x10 12 T°’ 5 

est 

K 18 /K 35 

0 . 1 14x10 19 T _1 ‘ ^exp(-7 1 , 600/T) 

50 

0 . 68exp(-7 1 , 600/T) 

0.57xl0 18 T _1 ’ 5 exp(-53,400/T) 

est 

K 22 K 69 

0.57xl0 18 T _1 " 5 exp(-33,800/T) 

est 

K 23 K 69 

0.57xl0 l8 T 1 ' 5 exp(-69,700/T) 

est 

0 .86exp(-69 , 700/T) 

0.57xlO J ' 8 T‘ 3 " 5 exp(-69,7O0/T) 

est 

K 72 

0.25x10 12 T°' 5 

est 

* 18*72 

0.25x10 12 T° ‘ 5 exp(-ll,400/T) 

est 

K 17 K 72 

0.30x10 15 T °* 5 exp(-16 , 1O0/T) 

est 

K 20 K 72 

0.30xl0 15 T’"°' 5 exp(-27,800/T) 

est 

CN] 

*4 

cr\ 

0-25x10 12 T° ’ 5 exp(-4600/T) 

est 

K 72 /K 16 

0.50x10 12 T°' 5 

est 

>0 

to 

AD 

0.50x10 L9 T -1 '° 

est 

K 9 K 72 

0.20x10 20 t' L '° 

est 

7 * 

Co 

O 



TABLE 3 (continued) 


82. N 

83. N 

84. C 

85. C 

86. N 

87. C 

88. N 

89. N 

90. C 

91. N 

92. 

93. 



Forward Rate Constant, 
cm /mole-sec 

Ref 

Equilibrium Constant 

Reaction 

N + M - N * + M. 

i. Z, ”*■ 

0.104xl0 l8 fc' 015 

es t 

K 5 K 69 

+ N + M 2 - N 2 * + M 2 

0.48x10 22 T _I ' 5 

es t 

K 82 

+ N + - CN* + Mj 

0.25xl0 15 

est 

K L1 K 86 

f N + M 2 - CN* + M 2 

0 .15xl0 16 

est 

K 85 

+ CN — CN* + N 

0.384xlO I6 T“°’ 5 exp(~36,900/T) 

est 

1.08esp(-36, 900/1) 

+ CN - CN* + C 

k 86 

est 

K 86 

+ CO -» CN* + 0 

0.57xl0 18 T _1 ’ 5 exp(-68 ) 800/T) 

est 

K 86 /K 19 

+ CO* - CN* + 0 

0.35xlO i5 T"° -5 exp(-900/T) 

est 

K 88 /K 7Z 

+ N 2 CN* +■ N 

0 . U4x 10 19 T _L * 5 exp(-53,800/T) 

est 

K 86^22 

+ C 2 - CN* + C 

0 . 10xl0 15 T~°' 5 exp(-ll , 300/T) 

est 

K 25 K 86 

+ CN* - N 2 * + C 

0 . 10x10 15 t“ 0 ’ 5 exp(-17 , 800/T) 

est 

K 69 /K 90 

+ NO - CN* + 0 

0.10xl0 1 5T‘°‘ 5 exp(-16 t 00/T) 

est 

K 24^ K 86 



possibilities are included in Table 3, and are written in exothermic form. 
Four possibilities were excluded because of spin conservation requirements 
or geometrical configuration. These are CO^ + C, CO^ + C -* 

CO + CO, CN + 0 2 - C0 2 + N, and C0 2 + N - CO + NO. 

11 % 3 

The basic estimated rate constant used here is 5x10 T cm /mole-sec, or 
about 1/6 the gas kinetic value. Many of the reactions considered can 
go to a number of different products, for example NO + C goes to either 
CO + N or CN + 0. The estimated rate constants for such reactions are 
chosen to be smalLer than the basic value by a factor equal to the number 
of possible sets of products. The reactions and rate constants used are 
given as 16 through 25 in Table 3. 

Two of the reaction rates listed above are an additional factor of 2 and 
5 slower, respectively, than the original estimate; a second is an addi- 
tional factor of 4 faster. This was done to improve the comparison 
between computed and measured values for non-equilibrium radiation (see 
Section 3.2.4). 

f. Ion Production and Loss . In air the rate of production of 

molecular ions behind strong shock waves is predicted satisfactorily 

by including only the associative ionization reactions, i.e., A + B “* 

+ 

AB + e + M. Because of the lack of knowledge concerning molecular 

-f- • 

impact ionization, i .e. , A^ + 1! -* A^ + e + M, this process is 
neglected here. Similarly electron impact ionization is neglected 
because of the uncertainty in electron temperature. If the analysis 

47 

carried out for monatomic gases is used, the electron impact rate is 
slower than the associative ionization rate for normal shock speeds up 
to 9 km/sec. Too little information is available on the electron tem- 
perature in an ionizing diatomic gas; hence the electron impact ioniza- 
tion process is neglected here. The rate constant is defined for reactions 
of the form A + A — A^ + + e by dlA^j/dt = -1/2 d[A]/dt = Aj . The 

reactions and rate constants used are given as 26 through 31 in Table 3. 
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g. Ion and Charge Transfer . In a three atom system there is a 
multitude of possible ion and charge transfer reactions involving atomic 
or molecular ions. Unfortunately none of these has been measured at low 
energies and we must rely on estimates. Only binary reactions involving 
one ion and one neutral are included here: reactions involving 2 ions 

are much less frequent because of the generally low degree of ionization. 


Most direct charge transfer reactions exhibit an appreciable energy 

48 

defect and, based on higher energy measurements, are expected to be 
much slower than gas kinetic. In a careful evaluation of the energy 
defect in direct charge transfer reactions in the C, N, 0 system, in- 
cluding low excited electronic levels, only two reactions were found 
which showed a sufficiently small energy defect (< 0.2 ev) to be con- 
sidered near-resonant. However ion charge is readily passed around 
through atom transfer reactions involving atomic or molecular ions for 
which there is no restriction on reaction rate other than the gas kinetic 

upper limit for exothermic reactions. The basic estimated rate constant 

11 ^ 3 

used for these reactions is 5x10 T cm /mole-sec, or about 1/6 the gas 
kinetic value, reduced by the number of possible sets of products. The 
reactions and rate constants used are given as 32 through 68 in Table 3. 

h. Electronic Excitation . The low electronic states (up to 2.5 ev) 
of all molecules and atoms are assumed to be in equilibrium with the 
ground states at the local translational temperature. This is justified 
on the basis that there are no known low excited states of or C0^ or 
A and that the low excited states of newly formed atoms and molecules 
are produced as rapidly as the ground states. Since some of these low 
excited states contribute appreciably to the heat capacity the required 
data is given in Table 4. 


The excitation mechanism of the electronic states is not well understood, 
but and CN impurity radiation is frequently observed immediately behmc 
the shock front in noble gas shock tube experiments, long before 
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TABLE 4 


Species 


CO 

N 


0 


2 


C 


2 


NO 


FACTOR BY WHICH THE INTERNAL ENERGY IS INCREASED 
DUE TO EXCITATION OF LOW ELECTRONIC STATES 

Factor 

1 

1 

§<-^“)exp(-ll,390/T) + |( 18 j,"° )exp(-18,990/T) 
I + | exp(-ll,390/T) + ~ exp(-18,990/T) 

( 2 - 7 ^ 5 ° )exp(-27,750/T) 

1 + 1 r exp(-27 , 750/T) 

1 


CO 

CN 


1 


2 (i3^30J)e x p(_i3 )3 oo/ T ) 

l - - 1- 2 exp(-13,300/T) 


i , 2. 5 ( 27, 700 /T) exp (-27 ,700/T) 
I + 2.5 exp(-27 ,700/T) 


. 555 ( 22^|50) exp (_22 f 850 /t) 
1 + I + .555 exp(-22 , 850/T) 


.555(14, 650/T)exp(- 14, 650/T) 
1 + .555 exp (-14, 650/T) 
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excitation or ionization of the noble gas is seen. This indicates that 

the rate of excitation of low electronic levels of molecules by neutral 

impact is very high. Large cross-sections have been observed for 

49 

quenching of upper electronic states by neutral impact, indicating high 
excitation rates for levels requiring only small electronic energy changes. 
This is extrapolated to low electronic levels of molecules, which require 
no change m total spin, by recognizing that the energy defect can be sup- 
plied through vibrational excitation of the molecule. Atoms in low 
excited states are then produced by atom transfer reactions. Low elec- 
tronic states of molecules which require a change in electronic spin may 
be excited primarily by free atom impact exchange reactions, as found. for 
the high electronic levels of Because most of the low electronic 

levels significant in radiation heating by the CNO system involve no spin 
change from the ground state, it was originally planned that the quasi- 
equilibnum assumption would be used. However, it was found that the 

results of this assumption did not fit the experimental data on CN violet 

2 +■ 

emission. Subsequently a finite rate of production of the B £ state of 
CN was introduced into the reaction rate system. 

In addition, the upper levels of and CO, which require a change m 

spin from those of the ground states, cannot be assumed to be in 

equilibrium with the ground states. The only direct experimental 

evidence available^ is for the reaction N + N~ -* N * + N for which 
18 -3/9 3 z / 

k = 1.14x10 T ~ exp(-71, 600/T) cm /mole-sec. This rate constant was 

found to be at least 100 times higher than that for the reaction 

N 2 * ^2 ^2* + N 2 " n0t k nown other free atoms, 0 for example, 

will cause electronic excitation of even though a direct atom 
transfer does not occur, but it is assumed here that they will not. 

In making estimates, all atom transfer reactions which can produce 
excited N or CO without violating the spin conservation rule are 
assumed to be at least as effective as N is in exciting N^. The ex- 
ponential is changed to the actual energy defect; the pre-exponential 


- 57 - 



is assumed to be the same as for N in exciting N^ if the energy defect 
exceeds 3 ev, a value 10 times larger is used if the energy defect is 
less than 3 ev. Excitation by neutral molecules is not included because 
of the lack of direct experimental evidence. Excitation by electron 
impact, although very efficient, cannot readily be included because of 
the lack of information concerning the electron temperature. Excited 
states above the first of ^ and CO are assumed to be in equilibrium 
with the first, denoted by N^* and CO*, at the local translational tem- 
perature. Excitation transfer between N^* and CO* is also included, 
assuming 1/4 the gas kinetic rate. 

Originally the rate constants for the atom transfer reactions that 
produce CN* were estimated to be the same as for N in exciting N 2 * 

However with these estimates, we found that the peak m CN* occurred 
much later than is seen experimentally. To obtain a good fit between 
theory and experiment, the reactions producing CN* were assumed to have 
rates approaching the gas kinetic value. 

Dissociation of excited molecules is also an important process, particularly 
near peak radiation, and is included in the reaction rate system for N^*, 
CO*, and CN*. The reactions and rate constants used are given as 69 
through 93 m Table 3. 

i. Equilibrium Constants . Not all of the reactions listed above 
are independent. A set of independent reactions which uniquely defines 
the equilibrium concentration of all species is formed. Equilibrium 
constants for this set are obtained by fitting the analytic expression 
aT b exp(c/T) to the results of the Aeronutronic equilibrium species com- 
puter program (Section 2,3.1). These expressions reproduce the computer 
results within a few percent over the temperature range 4000°K £ T 
10,000°K. Equilibrium constants for the remaining reactions are com- 
puted from combinations of the basic set. The values used m the 
computer program are listed in Table 3. 
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3.2.2 COLLISION LIMITING 


Some of the excited electronic states responsible for radiative heating 
in planetary entry may be de-excited by radiation more rapidly than 
they are excited by collision. This is known as "collision limiting". 
It requires that the radiative loss be included as a term in the reac- 
tion rate equation for the excited states affected by the radiation. 
Thus 


dN* 

dt 




Vj 



AN* 


(32) 


if the radiation is optically thin. 


If the radiation is optically thick, the radiation loss term is reduced 
by absorption. However this leads to a difficult integration problem 
because the solution at any point depends on the final solution at all 
points. To avoid this difficulty the radiation loss for any optically 
thick transition often can be neglected. As will be seen this difficulty 
does not occur for the present calculations. 

In order to determine which transitions are subject to collision limiting, 

the ratio AN */E Q Y.y is computed for the specific conditions of 
ij P ij 1 J 

interest. If it approaches or exceeds L, collision limiting must be 

included. For the present case and equilibrium conditions at the 

stagnation point, this ratio is about 1 for CN violet, 0.01 for CN red, 

0.05 for C0(4+), and 0.001 for CO (vib-rot) . This ratio is about the same 

at conditions of peak non-equilibrium radiation. From this, we conclude 

that collision limiting is significant only for CN violet emission. Since 

the CN violet emission is essentially optically thin for the present cases 

the radiation loss term can be included directly in the rate equation for 

7 -1 

CN*. This was done for the present calculations, using A = 1.06x10 sec 
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3.2.3 INTEGRATION SCHEME 

A Runge-Kutta numerical integration scheme, modified by Treanor,^ is 

used to compute the non-equilibrium species concentrations. Most of the 

details of the numerical method are identical to those described by 

Treanor"^ and are not repeated here. However the method by which the 

step size is chosen is considerably different, offering an appreciable 

gain in computer time. The allowable step size at each point is taken 

as a definite fraction of the shortest of the relaxation distances of 

each species. Based on experience a value for this fraction of 0.5 or 

less provides sufficient accuracy. However the step size must be reduced 

when one or more species approaches equilibrium or quasi-equilibrium con- 

52 

ditions. Under these conditions the rate equations become "stiff" and 
the modified Runge-Kutta numerical integration scheme leads to divergent 
oscillations if the step size is too large. Tests are applied which re- 
quire that the species production rate computed within the modified Runge- 
Kutta sequence must not exceed a value which would result m a relaxation 
distance shorter by a prescribed factor than the relaxation distance used 
m selecting the step size. If this test is failed, the step size is re- 
duced by a factor of 2. If this is passed on 2 succeeding steps the step 
size is increased by a factor of 2 for the next step, and on up to the 
original value of 0.5 times the minimum relaxation distance. 

The method permits integration of the non-equilibrium equations to within 
a few percent of equilibrium, requiring about 10 minutes on the Philco 
2000 computer for the CO^ system containing 13 species and 31 reactions, 
and 20 to 30 minutes for the N 2 ~C 02 “A system containing 24 species and 
93 reactions . 

3.2.4 COMPARISON WITH EXPERIMENT 

The reaction rate system described in the preceding paragraphs has been 
evaluated by calculating non-equilibrium radiation profiles behind normal 
shocks in CX^-^-A mixtures and comparing the results with experimental 
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measurements. In Section 3 .4 a description is given of the methods used 
in computing non-equilibrxum radiation profiles from the species profiles 
that are obtained from integrating the reaction rate equations. 

Four comparisons are shown m Figures 9, 10, 11, and 12. It is seen 
that the agreement is generally good. The CN violet profile behind a 
20,400 fps shock in a mixture of 25% CO^, 75 % N^ at 0.175 torr shows 
particularly good agreement (Figure 9). Since these conditions are 
close to those of the present flight case, very little extrapolation 
error is introduced. 

The CN red profile shown in Figure 10 is also in satisfactory agreement 
with the experimental data. The theoretical curve shows a somewhat 
sharper overshoot than the data shows: this is probably due m part 

to the finite slit width and the shock curvature present in the ex- 
periment. No other CN red comparisons are available. 

The CN red and CN violet profiles of Figures 10 and 11 indicate that the 
B £ excited state of CN is not in thermal equilibrium with the ground 
state. The peak experimental ratio of I to I is 4 for CN red and 5 for 
CN violet at t = 0.4 psec. If the excited states are both in equilibrium 
with the ground state, the translational temperature must be within 6% 
of its equilibrium value to yield I/I ratios this close to each other. 
However at the experimental peak, the translational temperature is pre- 
dicted to be 28% higher than the equilibrium value. For comparison, the 
CN violet emission predicted on the basis of equilibrium with the ground 
state is shown m Figure 11. The non-equilibrium profile fits the data 
more closely. 

The CN violet comparison shown in Figures 9, 11 and 12 are satisfactory 
in terms of computing the integrated non-equilibrium radiation. However 
the differences in profile shapes are sufficient to indicate that the 
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FIGURE 9. COMPARISON OF COMPUTED NON-EQUILIBRIUM CN VIOLET 

RADlAXION p PROFILE WITH THAT OBTAINED EXPERIMENTALLY 
BY ARNOLD 53 BEHIND A NORMAL SHOCK TRAVELING AT 20,400 FPS 
THROUGH A 25% C0 2 , 75% N 2 MIXTURE AT AU INITIAL PRESSURE 
OF 175 MUG. 




1-0 2.0 
t, n sec 


FIGURE 10. COMPARISON OF COMPUTED NON-EQUILIBRIUM CN RED RADIATION 
PROFILE WITH THAT OBTAINED EXPERIMENTALLY BY THOMAS AND 
MENARD55 BEHIND A NORMAL SHOCK TRAVELING AT 24,600 FPS 
IN A 9 % C02, 00% N2, 1% A MIXTURE AT AN INITIAL PRESSURE 
OF 250 jUHg. 





t , jU sec 


FIGURE 11. COMPARISON OF COMPUTED NON-EQUILIBRIUM CN VIOLET 

RADIATION PROFILE WITH THAT OBTAINED EXPERIMENTALLY 55 
FOR THE CONDITIONS OF FIGURE 30. DASHED CURVE IS 
THE NON -EQUILIBRIUM CN VIOLET RADIATION PROFILE 
COMPUTED ASSUMING THE RADIATING STATE IS IN 
EQUILIBRIUM WITH THE GROUND STATE OF CN AT THE 
NON-EQUILIBRIUM TRANSLATIONAL TEMPERATURE 







reaction rate system used here is not fully satisfactory. Significant 
improvements can be made in the future only with a great deal more 
experimental and analytical study. 

The approach and the reaction rate constants adopted here differ 

53 

significantly from those used by McKenzie. In his work the excited 
states of CN are assumed to be in equilibrium with the ground state. 

This makes it difficult to fit both CN red and CN violet data. His 
comparison with CN violet data is excellent, but no comparison with CN 
red is made. His prediction of the CN red peak is a factor of about 2.5 
less than that using our reaction rate system. This question bears 
further study. 

McKenzie's choice for the rate constant of one of the critical reactions, 
CO + N -* CN + 0, is less than a factor of 2 smaller than ours, showing 
good agreement. However his choice for + C — CN + N is a factor of 
7 smaller than ours at 9000°K. In addition his choice for the raLe of 
dissociation of CN exceeds our estimate by a factor of 15 (at 9000°K) 
for 0, C, „N } NO, and CN as collision partner and by a factor of 90 
(at 9000°K) for N^: CO, 0^, and as collision partner. McKenzie's 

choice appears quentionable until further support is obtained through 
direct experimental measurement. 

In summary, although the present reaction rate system is not wholly 
satisfactory, it can be expected to provide adequate accuracy for flight 
calculations of non-equilibrium radiative heat transfer, 

3.3 NON- EQUILIBRIUM CONVECTIVE HEAT TRANSFER ANALYSIS 

The convective heat tiansfer analysis utilized for non-equilibnum flow 
calculations is the same as described m Section 2.2 for equilibrium 
flows. The assumption made is that the heat transfer for frozen, non- 
equilibrium and equilibrium chemistry boundary lasers is identical even 
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for Lewis numbers slightly different than unity. This assumption is 
valid if (1) the edge of the boundary layer flow properties are at 
equilibrium, (2) the boundary layer gases are dissociated species that 
can be approximated as a binary mixture and (3) the wall is fully 
catalytic. The assumption of equilibrium edge properties is valid 
for the sphere cap, since the non-equilibrium property layer extends 
from the bow shock to a distance of only about 30% of the shock layer 
thickness ; t thus the properties at the boundary layer edge are very close 
to equilibrium. In the sphere-cone case the non-equilibrium region 
extends throughout the shock layer, but the boundary layer edge properties 
are not far from equilibrium conditions. 

The dissociated gas in the boundary layer is accurately approximated by 

a binary mixture. Because of the binary nature of the boundary layer gas, 
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the convective heat transfer is independent of the chemical kinetics. 

Radiation cooling of the inviscid gas is insignificant in the sphere-cone 
case, but must be included m the sphere-cap calculation. The radiation 
cooling effect is included m the stagnation-point convective heat transfer 
calculation by using the actual radiation-cooled boundary layer edge 
enthalpy rather than the adiabatic stagnation enthalpy. The convective 
heat transfer distribution over the whole body is thereby reduced by this 
correction factor. This is consistent with the result that radiation 
cooling reduces the boundary layer edge enthalpy by a factor that is 
nearly the same over the whole sphere-cap body. 

3.4 NON-EQUILIBRIUM RADIATIVE HEAT TRANSFER ANALYSIS 

3.4.1 NON- EQUILIBRIUM RADIATION 

The non-equilibrium species and excited state populations obtained as 
outlined in Section 3.2 are combined with the equilibrium radiation 
program described in Section 2.3 to provide non-equilibrium volume 
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radiation emission. The chemical model used in computing the excited 

state populations has been described in Sections 3.2. 1 and 3.2.2. The 

upper excited states of and CO are assumed to be in equilibrium with 

the lowest excited state, N *, and CO*, at the local non-equilibrium 

1 2 

translational temperature. The A il state of CN is assumed to be in 

2 + 

equilibrium with the ground state. The B £ state is computed directly 
as CN* m the integration of the reaction rate equations. The non- 
equilibrium excited state populations and the non-equilibrium temperature 
are used as inputs to the equilibrium radiation program to compute spectral 
absorption coefficients and non-equilibrium radiation emission. In this 
model the populations of the vibrational and rotational states of the 
excited band radiators are assumed to be in local thermodynamic equilibrium 
with the non-equilibrium translational temperature. This assumption is 
justified on the basis of shock tube measurements which show that the 
vibrational and rotational degrees of freedom are equilibrated immediately 
behind normal shockwaves for shock speeds above 13,000 to 15,000 ft/sec. 
Free-bound emission and absorption is computed on the basis of non- 
equilibnum species compositions and excited state populations in 
equilibrium with the nearest continuum of states at the local trans- 
lational temperature. The infrared vibration-rotation band emission of 
CO and NO and CN is assumed to be at the local non-equilibrium trans- 
lational temperature; this emission is only weakly dependent on 
temperature . 

Comparisons between calculations using the Aeronutronic non-equilibrium 
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radiation program and the shock tube measurements of Thomas and Menard 
53 

-and Arnold are shown m Figures 9, 10, 11, and 12 and are described in 
Section 3.2.4. These m tensity-time profiles indicate that the overall 
behavior of non-equilibrium radiation from CN is satisfactorily reproduced 
by the computer program. No comparison such as this can be made for the 
C0(4+-) band, because there is no published data. For this reason, the 


- 68 - 



C0(4+) band was assumed to reach its equilibrium value rapidly, buL to 
exhibit no overshoot. This assumption may be in error by one to two 
orders of magnitude. 

When the non-equilibrium integration along streamlines approaches 

equilibrium it is economical to shift from non-equilibrium radiation to 

equilibrium radiation. This is done at a point on the streamline for 

which the non-equilibrium radiation reaches within 10% of the local equi- 
$ 

librium value. When this shift is made on each streamline it is retained 
throughout the remainder of the flow field. Subsequent expansion of the 
streamline is treated as equilibrium, an approximation that is well 
justified in the present calculations on the basis of the non-equilibrium 
analysis . 

3.4.2 SELF- ABSORPTION OF NON- EQUILIBRIUM RADIATION 

The method used to account for self-absorption of non-equilibrium radiation 
is the same as that used for equilibrium radiation (see Section 2.3.3). 

This method involves the integration of the radiation emission out to a 
distance of one absorption length from the surface point. This is done 
m each direction along which the integration is made. The non-equilibrium 
absorption coefficient is the same as the equilibrium value at the same tem- 
peration and pressure, since the absorption arises from the ground electronic 
state of each molecule, and we have assumed equilibration of the rotational 
and vibrational states with the translational temperature. Fortunately 
the self- absorption corrections are small, affecting the CO fourth positive 
band strongly and the CN violet band slightly. The CO band contributes 
less than 107, of the radiation m the present non-equilibrium cases. 

3.4.3 RADIATION- INV IS Cl D FLOW COUPLING 

The cooling of the inviscid flow due to radiation Losses is included m 
these calculations by the following approximate method. Throughout most 
of the non-equilibrium region the effect reduces the enthalpy very little 
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and affects the total emission from the non-equilibrium layer by no more 
than 107<>. Because of the difficulty of including this effect in the 
integration of the rate equations, it has been neglected here. However 
beyond the point on each streamline at which the radiation is assumed 
to be the equilibrium value the effects of radiation cooling are included. 
This is done by utilizing the values of pressure, velocity, and density 
obtained from the flow field calculation (which neglected radiation 
cooling) and computing a first order correction to the enthalpy and the 
volume emission. This correction is found by integrating the (corrected) 
radiation loss along streamlines. In this integration self-absorption is 
included, although only minor effects arise m the present cases. 
Corrections to the volume emission are as large as 307„ m the sphere- 
cap case. Numerical results are given in Section 4.2. 

3.4.4 RADIATION-VISCOUS LAYER COUPLING 

This is included m the same manner as described for equilibrium radiation 
m Section 2.3.5. A portion of the C0(4+) band is self-absorbed within 
the thermal boundary layer thickness for the sphere-cap case and is there- 
fore omitted m the radiation integration. Radiation absorption m the 
boundary layer has negligible effect on the convective heat transfer. 

3.4.5 GEOMETRIC INTEGRATION 

This is carried out m much the same manner as in the equilibrium case 
described in Section 2.3 6. However the integration over the non- 
equilibrium layer for the sphere-cap case was simplified by determining 
the total emission from the non-equilibrium slab at each radial position. 
This was then used at each intersection of the viewing rays with the non- 
equilibrium slab when carrying out the geometric integration. 
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SECTION 4 


NUMERICAL RESULTS 


The flow field and the convective and radiative heat transfer have been 
computed for equilibrium flow over a sphere-cone and for non-equilibrium 
flow over a sphere-cone and sphere-cap for the following flight condition 


Mach Number 
Flight Velocity 
Ambient Density 
Ambient Temperature 
Ambient Composition 
Angle of Attack 
Base Diameter 


31 

19,600 ft/sec 

__ C. O 

10 slugs/ft 
180°R 

70% N , 30% C0 2 
0 ° 

12 ft 


Results of the equilibrium calculations are presented m Section 4.1; 
results of the non-equilibrium calculations are presented in Section 4 2. 


4.1 EQUILIBRIUM FLOW OVER A SPHERE-CONE 


4.1.1 EQUILIBRIUM SPHERE-CONE INVISCID FLOW PROPERTIES 
It was originally antic] pated that the mviscid flow field calculations 
would be made with a first approximation integral relation computer 
program. In order to verify the accuracy of the integral method, the 
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theoretical shock shape for a 60° half-angle cone m a Mach 9, Y = 1.4, 
air flow was computed and compared with experiment. The results are 
shown in Figure 5. It is seen that the first approximation solution 
provides an accurate calculation of the shock shape. It was determined, 
however, midway through the present program that real gas effects (y- 1«1) 
significantly influence the character of the 60° cone flow field, and 
that a higher order solution of the inviscid conservation equations is 
required to obtain an accurate shock shape and shock layer properties. 

It was decided that the Philco-Ford time dependent finite difference 
program should be used to make the inviscid flow field calculations. 

The sphere-cone shock shape obtained from both the time dependent finite 
difference equilibrium solution and a low Y (Y = 1.18) first approximation 
integral relation solution are shown m Figure 13. Although the integral 
solution is accurate in the sphere cap region, the shock layer thickness 
in the conical flow region is overestimated by as much as a factor of 2 
in some places. The first approximation integral relation solution is 
accurate for blunt body flow field calculations in which the shock layer 
flow is predominantly subsonic. When a large portion of the shock layer 
is supersonic as it is m the case of the equilibrium blunt 60° cone 
problem, however, the first approximation exaggerates the influence of 
the body properties on the shock layer calculation and overestimates the 
shock detachment distance. 

The time dependent solution indicates that the equilibrium shock layer 
in the conical region is predominantly supersonic and that the com- 
pression waves at the sphere cone junction travel without significant 
dispersion outward to the bow shock. The interaction produces a con- 
cave region m the bow shock where it regains strength. This concavity 
is observed in the shock shape plot presented in Figure 13. The sub- 
sequent reflections of the compression wave are negligible, but 
theoretically there is an entire series of alternating compression 
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TABLE 5 

EQUILIBRIUM INVISCID SHOCK LAYER PROPERTIES ON SPHERE- CONE 


Body 

Station 

Arc Length 
on Surface, 
Ft 

Standoff 

Distance, 

Ft 

Ft 

Ft 

Pressure 

Psia 

,Enthalpy 
Kcal/100 GM 

Spaed, 

10 ft/sec 

Radiatiog , 
watts/cm ster 

1 .Shock 

.0785 

.0677 

-.0622 

.0825 

2.33 

338 

.211 

.920 

Midpoint 



-.0282 

.0804 

2.44 

341 

.111 

.990 

Body 



+.00513 

.0783 

2.53 

342 

.0497 

1.038 

2 

.235 

.0694 

- .0422 

.247 

2.35 

332 

.313 

.878 




- .00786 

.240 

2.40 

335 

.255 

.930 




.0256 

.234 

2.45 

339 

.183 

.980 

3 

.392 

.0733 

- .00297 

.408 

2.07 

301 

.612 

.538 




.0320 

.396 

2.17 

321 

.442 

.723 




.0661 

.385 

2.20 

334 

.28-7 

.828 

4 

.549 

.0733 

.0605 

.562 

1.91 

269 

.817 

.306 




.0937 

.5 '.6 

1.92 

299 

.628 

.480 




.126 

.5_ J 

1.95 

326 

.388 

.673 

5 

.706 

.0767 

.134 

.706 

1.85 

251 

.911 

.246 




.167 

.687 

1.84 

282 

.738 

.362 




.200 

.668 

1.83 

315 

.496 

.555 

6 

.863 

.0821 

.207 

.845 

1.87 

247 

.931 

.204 




.243 

.824 

1.90 

273 

.793 

.332 




.279 

.804 

1.93 

309 

.550 

.546 

7 

1.021 

.0889 

.280 

.985 

1.91 

248 

.922 

.216 




.319 

.962 

1.94 

267 

.825 

.304 




.357 

.940 

1.97 

302 

.607 

.508 

8 

1.346 

.1083 

.426 

1.277 

1.95 

255 

.889 

.254 




.473 

1.249 

1.98 

261 

.857 

.287 




.520 

1.222 

1 99 

284 

.726 

.410 



Table 5 (continued) 


i 

->j 

■c- 

i 


Arc Length 
Body on Surface, 
S tation Ft 

9 1.841 

10 2.336 

11 2.830 

12 3.325 

13 3.820 

14 4.314 

15 4.809 

16 5.304 


Standoff 
Distance , 
Ft 

.1412 

.1805 

.2127 

.2536 

.2895 

.3251 

.3671 

.3979 

.441 


X, 

Ft 

y, 

. Ft 

Pressure, 

Psia 

Enthalpy Speed, 

Kcal/100 GM 10 ft/sec 

Radiatio^, 
watts/cm ster 

.644 

1.722 

1.96 

258 

.876 

.262 

.706 

1.686 

1.99 

260 

.863 

.284 

.767 

1.651 

2.00 

275 

.779 

.362 

.860 

2.168 

1.97 

2 60 

.862 

.282 

.938 

2 .123 

2.00 

2 64 

.862 

.308 

1.015 

2.079 

2.00 

270 

.808 

.340 

1.076 

2.615 

1.97 

259 

.868 

.277 

1.170 

2.561 

1.99 

260 

.864 

.285 

1.262 

2.507 

2.00 

267 

.825 

.323 

1.293 

3.061 

1.97 

260 

.861 

.282 

1.402 

2.998 

1.99 

260 

.863 

.285 

1.509 

2.939 

2 .00 

265 

.836 

.313 

1.508 

3.508 

1.97 

' 260 

.863 

.282 

1.634 

3-435 

1.99 

260 

.864 

.285 

1.757 

3.364 

2.00 

264 

.844 

.308 

1.726 

3.953 

1.96 

250 

.871 

.275 

1.866 

3.872 

1.99 

259 

.866 

.280 

2.004 

3.793 

2.00 

262 

.850 

.298 

1.942 

4.400 

1.97 

261 

.858 

.288 

2.098 

4.310 

1.99 

260 

.865 

.285 

2.251 

4.221 

2.00 

261 

.856 

.292 

2.159 

4.845 

1.96 

257 

.876 

.267 

2.331 

4.746 

1.99 

259 

,869 

.280 

2.499 

4.649 

2.00 

260 

.860 

2.37 

2.379 

5.289 

1.97 

259 

.866 

.277 

2.567 

5.181 

1.99 

259 

.868 

.280 

2.746 

5.078 

2.00 

260 

.864 

.287 


17 


5.798 



Table 5 (continued) 


I 

Ln 

l 


Body 
S tatxon 

18 


19 


20 

21 


22 


23 

24 


25 


Arc Length 
on Surface, 
Ft 

S tandof f 
Distance , 
Ft 

Ft , 

y. 

Ft 

6.124 

.468 

2.522 

5.584 



2.721 

5.468 



2.909 

5.360 

6.281 

.478 

2.590 

5.726 



2.798 

5.606 



2.938 

5.496 

6.438 

.487 

2.657 

5.868 



2.881 

5.739 



3.066 

5.632 

6.595 

.520 

2.759 

6.069 



2.982 

5.895 



3.153 

5.762 

6.752 

.608 

2.909 

6.348 



3.113 

6.073 



3.263 

5.872 

6.909 

.784 

3.117 

6.663 



3.282 

6.241 



3.397 

5.9 50 

7.065 

1.072 

3.434 

7.052 



3.502 

6.400 



3.547 

5.990 

7.222 . 

1.617 

3.945 

7.532 



3.787 

6.545 



3.703 

5.990 

7.378 

2.581 

4.788 

8.054 



4.157 

6.646 



3.853 

5.950 


Pressure, 

Psia 

Enthalpy 
Kcal/100 GM 

SpjiGdj 
10 ft/sec 

Radiatiop 
watts/ era . 

1.98 

261 

.860 

.289 

2.00 

259 

.868 

.282 

2.00 

250 

.866 

.282 

1.97 

261 

.860 

.288 

1.99 

259 

.868 

.280 

2.00 

259 

.867 

.282 

1.97 

261 

.857 

.288 

1.99 

259 

.870 

.280 

2.00 

259 

.869 

.282 

1.93 

257 

.879 

.258 

1.88 

256 

.883 

.244 

1.70 

253 

.899 

.200 

1.87 

250 

.913 

.215 

1.62 

250 

.911 

.177 

1.16 

245 

.940 

.102 

1.66 

233 

.995 

.123 

1.29 

244 

.944 

.115 

. .713 

238 

.971 

.047 

1.42 

199 

1.137 

.0364 

.970 

237 

.975 

.0680 

.435 

239 

.967 

.0262 

1.17 

186 

1.190 

.0171 

.690 

230 

1.007 

.0366 

.264 

23'8 

.973 

.0135 

.821 

193 

1.162 

.0151 

.449 

.157 

221 

234 

1.047 

.991 

.0171 

.0064 


26 



RADIAL COORDINATE 
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and expansion waves along the layer. In the higher y cases, the subsonic 
shock layer disperses the compression and there is no concavity to the 
shock. At higher y's, the concavity has been observed on sharper cones. 

The flow over the 60° sphere-cone in air (y = 1.4) shown m Figure 5 con- 
tains a subsonic shock layer with a sonic line running approximately normal 
to the cone surface and starting at the outer corner of the cone. In con- 
trast, the equilibrium flow of a 307„ CC^, 70% mixture over a similar 
body shape displays a shock layer that is almost entirely supersonic. 

The velocity on the entire conical surface is subsonic, but, as seen in 
Figure 15, it rapidly becomes supersonic away from the surface. The sonic 
line runs from the rounded shoulder of the base roughly parallel to the 
cone surface at about 10% of the standoff distance, finally crossing the 
shock layer below the sphere-cone junction to reach the bow shock. 

In Figure 14, the cone surface pressure as a function of the ratio of body 
surface distance to base radius is shown. The two prominent features of 
the pressure distribution a re the rapid decrease in pressure on the 
spherical po_rtion of the sphere-cone and the pressure recovery at the 
sphere-cone junction to a pressure plateau in the conical region. The 
pressure distribution in the nose region is similar to a hemispherical 
shock layer result and the pressure plateau m the conical region is 
similar to a pointed cone supersonic flow result. In Figure 14 it is seen 
that the Newtonian flow model, which is a low y limit of the hypersonic 
inviscid flow field equations, provides a useful approximation to the 
detailed flow field pressure distribution result. 

In Table 5, the properties of the shock layer gas on 26 normals to the 
body surface are tabulated. For reference, the location of the body 
normals, in which the shock layer properties are specified, are shown on 
the body in Figure 13- The arc length, s, along the body to the inter- 
section with each normal and the shock standoff distance on each normal 
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RESSURE 



FIGURE 14. EQUILIBRIUM SURFACE PRESSURE ON 



SPHERE-CONE 




FIGURE 15. FLUID VELOCITY PROFILES ON BODY NORMALS OF 
THE SPHERE- CONE 
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are given in the second and thud columns. The axial and radial 
coordinates of the intersections of each normal with the shock, with 
the half-way point, and with the body surface are specified m the 
vertical grouping in the next two columns. The origin of the coordinate 
system is at the nose of the sphere. The pressure, enthalpy, temperature, 
fluid velocity and radiation flux at the three points on each normal are 
listed in the subsequent columns. 

In this tabulation, the presence of an actual inviscid entropy layer on 
the conical portion has been neglected. This layer contains the fluid 
that entered the shock, layer in the blunt nose region An estimate of 
the width of the layer can be obtained from the fine mesh computation 
discussed earlier for y = 1.1, M = 20. A typical velocity 
profile along a normaL to the conical portion is presented in Figure 15. 
The normal, is located near the onset of the region of constant pressure 
on the body surface, where the entropy Layer has its largest extent. It 
is seen that the layer comprises at most 257. of the entire shock layer. 
Since the viscous boundary layer can be expected to comprise 157, of the 
shock layer, width, there should be only a minor entropy layer. The 
coarse mesh, which averages quantities along each one-third of the normal, 
adequately portrays the flow encountered by the outer edge of the boundary 
layer . 


4.1.2 EQUILIBRIUM SPHERE-CONE CONVECTIVE HEAT TRANSFER 
The equilibrium sphere-cone laminar convective heat flux versus distance 
along the body is shown in Figure 16. For comparison purposes the heat 
flux obtained from modified Newtonian pressure distribution and normal 
shock entropy is also plotted. In the spherical region of the body, 
the equilibrium laminar convective heat transfer exceeds the modified 
Newtonian result because of the vorticity interaction effects. In the 
region aft of the sphere-cone junction, the detailed flow field solution 
predicts higher heat transfer because the boundary layer edge velocity 
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HEATING RATE, BTU/FT -SEC 



JUNCTION DISTANCE ALONG THE' BODY, S/R_ 


FIGURE 16. SPHERE-CONE EQUILIBRIUM CHEMISTRY CONVECTIVE 
AND RADIATIVE HEAT TRANSFER DISTRIBUTIONS 
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conditions in this region are characterized by the conical shock entropy 
rather than the normal shock entropy. As a result the boundary layer edge, 
velocities and hence the heat transfer are higher than those given by the 
Newtonian normal shock entropy calculation. On the sphere, the combined 
effects of vorticity interaction (12%) and ambient composition (8%) 
increase the stagnation point heat transfer approximately 20%. over an 
equivalent high Reynolds number air calculation. 

4.1.3 EQUILIBRIUM SPHERE- CONE RADIATIVE HEAT TRANSFER 

The optically thin spectrally integrated radiative emission from volume 
elements of the shock layer gas is included m Table 5. The spectral 
distribution is given in Table 6 for the gas adjacent to the stagnation 
point and m Table 7 for the gas adjacent to point 5 on the sphere-cone 
(see Figure 13). The species concentrations for each of these points 
are also given in the respective tables. The major radiation arises from 
the CN violet and red bands and from the 4th positive band of CO. 

Minor contributions arise from the E-X and H-X bands of CN and the 
vibration-rotation bands of CO. 


Because self-absorption of the CN violet band may be greater than 
inferred from the course frequency interval used in obtaining the 
results of Table 6, a much finer frequency interval was also used. 
These results are given in Table 8, and show a higher maximum 
absorption coefficient by a factor of 2.5. However even with this 
value, there is negligible self-absorption of the CN violet band 
m the equilibrium sphere-cone flow field. 


The radiative heat transfer to the seven body points shown m Figure 
13 is given in Table 9. In the final column the predicted radiation 
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TABLE 6 


i 

00 
to 

1 


EQUILIBRIUM RADIATION FROM FLOW ADJACENT TO SPHERE-CONE STAGNATION POINT 

(T = 5735°K, P = 0.172 atm) 


Mole-Fraction for CN VLT = 0.37560-002 Mole-Fraction for CN RED = 0.37560-002 
Mole-Fraction for CN E-X = 0.37560-002 < Mole-Fraction for CN H-X = 0.37560-002 
Mole-Fraction for 00(4+) = 0.18280+000 


Band System Opacity, CM 

CN RED CN E-X CN H-X 


0.20000+004 

0.40000+004 

0.60000+004 

0.80000+004 

0.10000+005 

0.12000+005 

0.14000+005 

0.16000+005 

0.18000+005 

0.20000+005 

0.22000+005 

0.24000+005 

0.26000+005 

0.28000+005 

0.30000+005 

0.32000+005 

0.34000+005 

0.36000+005 

0.38000+005 

0.40000+005 

0.42000+005 

0.44000+005 

0.46000+005 

0.48000+005 

0.50000+005 

0.52000+005 

0.54000+005 


0. 

0. 

0 46148-014 
0.30425-012 
0.65861-011 
0.57831-010 
0.90560-010 
0 26066-008 
0.19210-006 
0.43475-005 
0.63090-004 
0.63 983-003 
0.62833-002 
0.13411-002 
0.35291-004 
0.31430-005 
0.56086-007 
0.93277-008 
0.11921-009 
0.38280-010 
0.20920-011 
0.30676-013 
0.16939-014 
0.24141-015 
0.11412-016 
0.67236-018 
0.43838-021 


0.14312-004 
0.93002-004 
0.28116-003 
0.40538-003 
0.55264-003 
0 30772-003 
0.16743-003 
0.65912-004 
0.20172-004 
0.54522-005 
0.13089-005 
0.27540-006 
0.46057-007 
0.66029-008 
0.78125-009 
0.72126-010 
0.48608-011 
0.21443-012 
0. 

0. 

0. 

0. 

0 . 

0. 

0. 

0. 

0 . 


0.49105-025 
0.74936-024 
0.91666-023 
0.10280-021 
0.11050-020 
0.11610-019 
0.12042-018 
0.12396-017 
0.12704-016 
0.12986-015 
0.13254-014 
0.L3616-013 
0.13876-012 
0.14142-011 
0.14422-01 0 
0.14714-009 
0.14984-008 
0.15266-007 
0.15585-006 
0.15861-005 
0.13578-004 
0.62911-004 
0.14935-003 
0.30665-003 
0.59595-003 
0.11861-002 
0 . 20201-002 


0.81808-028 
0.15086-026 
0.22371-025 
0.30441-024 
0.39716-023 
0.50658-022 
0.63792-021 
0.79735-020 
0.99223-019 
0. L2315-017 
0.15261-016 
0.18893-015 
0.23376-014 
0.28913-013 
0.35754-012 
0.44208-011 
0.54657-010 
0.67576-009 
0.83854-008 
0.10369-006 
0.12685-005 
0.10738-004 
0.40411-004 
0.86559-004 
0.17406-003 
0.36963-003 
0.94194-003 


0.72822-033 

0.22247-031 

0.55076-030 

0,12539-028 

0.27394-027 

0.58535-026 

0.12351-024 

0.25872-023 

0.53958-022 

0.11224-020 

0.23313-019 

0.48375-018 

0.10032-016 

0.20798-015 

0.43109-014 

0.89343-013 

0.18515-011 

0.38366-010 

0.79500-009 

0.16475-007 

0.34179-006 

0.70090-005 

0.12631-003 

0.13621-002 

0.52959-002 

0.12400-001 

0.24544-001 



Table 6 (continued) 


Wave Number 


Band System Opacity, 

cm" 1 


CM -i 

CN VLT 

CN RED 

CN E-X 

CN H-X 

CO (4+) 

0.56000+005 

0.13550-026 

0 . 

0.29352-002 

0.80197-003 

0.46001-001 

0.58000+005 

0.41882-032 

0 . 

0.41230-002 

0.11818-002 ' 

0.94080-001 

0.60000+005 

0.12946-037 

0 . 

0.48172-002 

0.15242-002 

0.16549+000 

0.62000+005 

0.40015-043 

0 . 

0.54684-002 

0.19009-002 

0.24254+000 

0.64000+005 

0.12369-048 

0 . 

0.38690-002 

0.22547-002 

0.34075+000 

0.66000+005 

0.38231-054 

0 . 

0.23734-002 

0.20347-002 

0.38584+000 


Wave Number, 
cm 

Total Absorption Coefficient, 
cm 

B(W,T), 2 
watts/cm -ster 

Emission,^ 
watts/cni - ster 

0.20000+004 

0.23336-003 

0.14620-001 

0.34119-005 

0.40000+004 

0.99729-004 

0.44110-001 

0.43990-005 

0.60000+004 

0.28117-003 

0.73381-001 

0.20632-004 

0.80000+004 

0.40538-003 

0.94660-001 

0.38374-004 

0 . 10000+005 

0.55264-003 

0.10548+000 

0.58292-004 

0.12000+005 

0.30772-003 

0.10663+000 

0.32812-004 

0.14000+005 

0.16743-003 

0.10047+000 

0.16821-004 

0.16000+005 

0.65915-004 

0.89711-001 

0.59133-005 

0.18000+005 

0.20364-004 

0.76780-001 

0.15636-005 

0.20000+005 

0.97998-005 

0.63492-001 

0.62221-006 

0.22000+005 

0.64399-004 

0.51032-001 

0.32864-005 

0.24000+005 

0.64011-003 

0.40051-001 

0.25637-004 

0.26000+005 

0.62834-002 

0.30801-001 

0.19353-003 

0.28000+005 

0.13411-002 

0.32379-001 

0.31220-004 

0.30000+005 

0.35292-004 

0.17329-001 

0.61159-006 

0.32000+005 

0.31432-005 

0.12731-001 

0.40016-007 

0.34000+005 

0.57626-007 

0.92445-002 

0.53272-009 

0.36000+005 

0.25291-007 

0.66436-002 

0.16803-009 

0.38000+005 

0.16514-006 

0.47306-002 

0.78119-009 

0.40000+005 

0. 17063-005 

0.33406-002 

0.56999-008 



Table 6 (Continued) 


Wave Number, 
cm 

Total Absorption_£oef flcient 9 
cm 

B(W, T) , 2 
watts/cm -ster 

Emission* 
watts/cm -ster 

0.42000+005 

0.15188-004 

0.23414-002 

0.35561-007 

0.44000+005 

0.80658-004 

0.16299-002 

0.13147-006 

0.46000+005 

0.31607-003 

0.11276-002 

0.35641-006 

0.48000+005 

0.17553-002 

0.77571-003 

0.13616-005 

0.50000+005 

0.60659-002 

0.53085-003 

0.32201-005 

0.52000+005 

0.13956-001 

0.36154-003 

0.50456-005 

0.54000+005 

0.27506-001 

0.24514-003 

0.67428-005 

0.56000+005 

0.49738-001 

0.16553-003 

0.82333-005 

0.58000+005 

0.99385-001 

0.11135-003 

0.11067-004 

0.60000+005 

0.17183+000 

0.74637-004 

0.12825-004 

0.62000+005 

0.24991+000 

0.49861-004 

0.12461-004 

0.64000+005 

0.34687+000 

0.33206-004 

0.11518-004 

0.66000+005 

0.39024+000 

0.22049-004 

0.86046-005 


Total Emission = 0.10376+001 watts/cm -ster 



TABLE 7 


EQUILIBRIUM RADIATION FROM FLOW ADJACENT TO BODY POINT 5 ON SPHERE-CONE 

(T*5342°K, P = 0.136 atm) 


Mole-Fraction for CN VLT = 0.16110-002 ' Mole-Fraction for CN RED*= 0.16110-002 

Mole-Fraction for CN E-X = 0.16110-002 Mole-Fraction for CN H-X = 0.16110-002 


Mole-Fraction for CO(4+) = 0.20980+000 

Wave Number, 

CM CN VLT CN RED 


Band System Opacities, CM 


.-1 


CN E-X 


CN H-X 


C0(4+) 


0.40000+004 
0.60000+004 
0.80000+004 
0.10000+005 
0.12 000+005 
0,14000+005 
0.16000+005 
0.18000+005 
0.20000+005 
0.22000+005 
0.24000+005 
0.26000+005 
0.28000+005 
0.30000+005 
0.32000+005 
0.34000+005 
0.36000+005 
0.38000+005 
0.40000+005 
0.42000+005 
0.44000+005 
0.46000+005 
0.48000+005 
0.50000+005 
0.52000+005 
0.54000+005 


0 . 

0.13256-014 
0.87773-013 
0.19116-011 
0.16889-010 
0.26356-010 
0.77793-009 
0.58768-007 
0. 13788-005 
0.21010-004 
0.22749-003 
0.23798-002 
0.48813-003 
0.12414-004 
0.10138-005 
0.17910-007 
0.30806-008 
0.37406-010 
0.13010-010 
0.69667-012 
0.10166-013 
0.60152-015 
0.88512-016 
0.43113-017 
0.26177-018 
0.10114-021 


0.32935-004 
0.10525-003 
0.15728-003 
0.21970-003 
0 1 L603-003 
0.61615-004 
0.23787-004 
0.70934-005 
0.18701-005 
0.44494-006 
0.93808-007 
0.15902-007 
0.23180-008 
0.28053-009 
0.26606-010 
0. 18454-011 
0.84607-013 
0.19589-014 
0. 

0. 

0 . 

0 . 

0 . 

0. 

0. 

0 . 


0.81258-026 
0.11653-024 
0.15368-023 
0.19472-022 
0.24158-021 
0.29627-020 
0.36097-019 
0.43816-018 
0.53070-017 
0.64198-016 
0.77604-015 
0.93769-014 
0.11327-012 
0.13682-011 
0 . 1662L-010 
0.20081-009 
0.24309-008 
0.29371-007 
0.35397-006 
0.34702-005 
0.17425-004 
0.43643-004 
0.94267-004 
0.19301-003 
0 .40350-003 
0.71894-003 


0.10980-028 

0.19363-027 

0.31431-026 

0.49034-025 

0.74917-024 

0.11315-022 

0.16980-021 

0.25386-020 

0.37872-019 

0.56430-018 

0.84020-017 

0.12505-015 

0.18606-014 

O.27681-0L3 

0.41178-012 

0.61254-011 

0.91117-010 

0.13588-008 

0.20216-007 

0.29609-006 

0.28210-005 

0.11376-004 

0.25640-004 

0.54188-004 

0.11940-003 

0.31936-003 


0.13753-033 
0.42056-032 
0.11865-030 
0.32197-029 
0.85602-028 
0.22504-026 
0.58786-025 
0.15301-023 
0.39740-022 
0.10309-020 
0.26723-019 
0.69245-018 
0.17938-016 
0.46464-015 
0.12034-013 
0.31167-012 
0 . 807 ] 6-011 
0.20904-00 9 
0.54134-008 
0.14020-006 
0.35745-005 
0.77036-004 
0.93021-003 
0.38329-002 
0.94365-002 
0.19582-001 



I 

cc 

*‘*sj 

I 



0.56000+005 

0. 

o. 

0 . 10803-002 

0.28494-003 

0.38632-001 

0.58000+005 

0. 

0. 

0.15562-002 

0-.431 17-003 

0.83060-001 

0.60000+005 

0. 

0. 

0.18525-002 

0.57078-003 

0.15327+000 

0.62000+005 

0. 

0. 

0.21381-002 

0.72359-003 

0.23249+000 

0.64000+005 

0. 

0. 

0.14899-002 

0.87851-003 

0.33524+000 

0.66000+005 

0. 

0. 

0.90220-003 

0.78894-003 

0.38462+000 

0.68000+005 

0. 

0. 

0.48915-003 

0.69110-003 

0.47296 +000 


Wave Number, Total Absorption Coefficient, B(W,T), ? Emission^ 


- I 

cm 

JL 

cm 

watts/cm" -ster 

watts/ cm -ster 

0.40000+004 

0.43746-003 

0.39350-001 

0.17214-004 

0.60000+004 

0.11704-003 

0.63781-003 

'0.74651-005 

0.80000+004 

0.15731-003 

0.79962-001 

0.12579-004 

0.10000+005 

0.21972-003 

0.86412-001 

0.18986-004 

0.12000+005 

0.11608-003 

0.84575-001 

0.98177-005 

0.14000+005 

0.61695-004 

0.77049-001 

0.47536-005 

0.16000+005 

0.23883-004 

0.66460-001 

0.15872-005 

0.18000+005 

0.72562-005 

0.54905-001 

0.39840-006 

0.20000+005 

0.33335-005 

0.43804-001 

0.14602-006. 

0.22000+005 

0.21500-004 

0.33956-001 

0.73007-006 

0.24000+005 

0,22759-003 

0.25695-001 

0.58480-005 

0.26000+005 

0.23799-002 

0.19051-001 

0.45339-004 

0.28000+005 

0.48828-003 

0.13879-001 

0.67770-005 

0.30000+005 

0.12617-004 

0.99589-002 

0.12566-006 

0.32000+005 

0.12398-005 

0.70518-002 

0.87427-008 

0.34000+005 

0.25483-006 

0.49352-002 

0.12577-008 

0.36000+005 

0.24834-006 

0.34184-002 

0.84892-009 

0.38000+005 

0.27713-006 

0.23459-002 

0.65010-009 

0.40000+005 

0.62899-006 

0.1596 6-002 

0.10042-008 

0.42000+005 

0.41577-005 

0.10785-002 

0.44840-008 



TABLE 7 (conintued) 


Wave Number, 
cm 

Total Absorption Coefficient, 
cm 

B ( W,T), 

watts/cm -ster 

Emission, 
watts/cnr -ster 

0.44000+005 

0.24073-004 

0.72355-003 

0.17418-007 

0.46000+005 

0.13232-003 

0.48244-003 

0.63834-007 

0.48000+005 

0.10504-002 

0.31985-003 

0.33597-006 

0.50000+005 

0.40803-002 

0.21095-003 

0.86076-006 

0.52000+005 

0.99597-002 

0.13847-003 

0.13791-005 

0.54000+005 

0.20621-001 

0.90484-004 

0.18658-005 

0.56000+005 

0.39998-001 

0.58886-004 

0.23553-005 

0 . 58000+005 

0.85048-001 

0.38176-004 

0.32467-005 

0.60000+005 

0.15569+000 

0.24661-004 

0.38395-005 

0.62000+005 

0.23535+000 

0.15878-004 

0.37368-005 

0.64000+005 

0.33760+000 

0.10191-004 

0.34405-005 

0.66000+005 

0.38632+000 

0.65217-005 

0.25194-005 

0.68000+005 

0.47414+000 

0.41621-005 

0.19734-005 


3 

Total Emission = 0.31484+000 watts/cm -ster 



TABLE 8 

DETAILED SPECTRAL DISTRIBUTION OF OPACITY 
OF CN VIOLET BAND 
(T = 5735°K, P = 0.172 atm) 

Mole-Fraction for CN VLT = 0.37560-002 


Wave Number, 
cm 


Opacity, 
cm 


0.21000+005 

0.21400+005 

0.21800+005 

0.22200+005 

0.22600+005 

0.23000+005 

0.23400+005 

0.23800+005 

0.24200+005 

0.24600+005 

0.25000+005 

0.25400+005 

0.25800+005 

0.26200+005 

0.26600+005 

0.27000+005 

0.27400+005 

0.27800+005 

0.28200+005 

0.28600+005 

0.29000+005 

0.29400+005 

0,29800+005 


0 . 

0 . 

0.70214-004 

0.17970-003 

0.61393-004 

0.49035-005 

0.38785-006 

0.12719-002 

0.16983-002 

0.21432-003 

0.16950-004 

0.13405-005 

0.16090-001 

0.14051-001 

0.11970-002 

0.94659-004 

0.74856-005 

0.13008-002 

0.49598-002 

0.39264-003 

0.31049-004' 

0.34936-005 

0.37815-004 
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TABLE 9 

EQUILIBRIUM SPHERE-COME RADIATIVE HEAT TRANSFER 

R^, watts/cm^ 

2 2 (with self-absorption) 

Position , watts/cm , watts/cm' 1 and boundary layer 

(Figure 13) (optically thin) ( with self-absorption) correction) 


1 


5.6 

5.4 

4 . 6 

2 


5.2 

5.1 

4.3 

3 


3.1 

3.0 

2.6 

4 


5.6 

5.2 

4.8 

5 


9.2 

8.1 

7.6 

6 


10.0 

8.5 

8.0 
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is given and includes the effects of self-absorption and viscous-inviscid 
coupling. It was found that only CO(4-t) emission was subject to significant 
self-absorption. The magnitudes of the self-absorption and boundary-layer 
coupling are seen by comparing the three columns. 

4.2 NON- EQUILIBRIUM FLOW OVER A SPHERE-CONE AND A SPHERE-CAP 

4.2.2 NON- EQUILIBRIUM INVISCID FLOW PROPERTIES 

a. Sphere Cone . The bow shock and shock layer streamline positions for 
the non-equilibrium sphere cone flow field are shown m Figure 17. For com- 
parison, the bow shock positions for the equilibrium solution and for a first 
approximation method of integral relation non-equilibrium solution are also 
shown. It is seen that the non-equilibrium shock layer detachment distance 
exceeds the equilibrium result by a factor of two in some regions. The dif- 
ference between these two cases arises from the lower shock layer densities 
and the more extensive region of subsonic flow in the non-equilibrium case. 

In the equilibrium calculation the temperature is lower throughout the flow 
field. This causes the gas density to be high and it also causes almost all 
of the flow to be supersonic at the sphere-cone junction. In the non-equilib- 
rium calculation, the temperature is higher and an extensive region of sub- 
sonic flow exists aft of the sphere-cone junction. Since the standoff distance 
varies inversely with density and velocity, the non-equilibrium standoff dis- 
tance is larger than the equilibrium value. Note also that the integral rela- 
tion solution for the non-equilibrium shock layer thickness somewhat exceeds 
the time-dependent result, but that the errors are not nearly as large as m 
the equilibrium case (Figure 13) . The reason that the integral solution is 
more accurate m the non-equilibrium case is that the sonic line in the one 
strip integral approximation, which extends linearly from the bow shock m 
the sphere-cone junction region to the rounded corner in the base region of 
the body, more closely approximates the actual flow conditions. 
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RADIAL COORDINATE 



IGURE 17. BLUNT SPHERE-CONE BODY, NON-EQUILIBRIUM SHOCK SHAPE, AND STREAMLINE 
POSITIONS 
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The body pressure distribution is shown in Figure 18. A Newtonian result 
is also shown for comparison. As in the case of the equilibrium flow, the 
two prominent features of the pressure distribution are the rapid decrease 
in pressure on the sphere-cone and the pressure recovery at the sphere- 
cone junction to a pressure plateau in the conical region. Even though 
the subsonic flow region extends beyond the sphere-cone junction, the 
non-equilibrium pressure distribution is very similar to the equilibrium 
result, although the actual pressure level is somewhat lower. Note also 
that the Newtonian model provides a useful approximation to the non- 
equilibrium pressure distribution. 

The flow properties throughout the inviscid region of the non-equilibrium 
sphere-cone flow field are given m further detail in Table 10. The 
notation follows that used in Table 5 and discussed in Section 4.1.1. 

Each body point is shown in Figure 17 . Data is given along rays normal 
to the body surface. The non-equilibrium species distributions are 
discussed in Section 4.2.3. 

b. Sphere-Cap . The sphere-cap bow shock and streamline field ob- 
tained from the time dependent finite difference solution are shown m 
Figure 19 . The bow shock is nearly concentric with the body except in 
the region of the corner. At the corner position, the sonic line extends 
m a normal direction from the body curving downstream slightly and then 
upstream to the bow shock. 

The inviscid result is based m a constant Y = 1.225, which corresponds 
to a shock layer with fully equilibrated rotational and vibrational energy 
modes and frozen free stream chemical species. In the course of performing 
the non-equilibrium streamline integrations, it was found that the non- 
equilibrium chemical species relaxation zone is only 307, of the shock 
layer thickness. Thus more than two-thirds o£ the shock layer is nearly 
m thermochemical equilibrium and the assumption of frozen chemical species 
shock layer is in error. This error in shock layer chemistry model causes 
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PRESSURE, 



FIGURE 18. NON-EQUILIBRIUM SURFACE PRESSURE ON THE SPHERE-CONE 
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TABLE 10 

NON-EQUILIBRIUM INVISCID SHOCK LAYER PROPERTIES ON THE SPHERE-CONE 


Body 
S tation 

Arc Length 
on Surface, 
Ft 

Standoff 

Distance, 

Ft 

' X, 

Ft 

1 .Shock 
.625 e 
.375 e 
Body 

.0785 

0.107 

-.1015 

- .0605 

- .0338 
.0051 

2. 

.235 

0.108 

-.0799 

.0393 

.0130 

.0256 

3. 

.392 

0.110 

- .0365 
.1835 
.0280 
.06612 

4. 

.549 

0.110 

.02227 

.06221 

.0895 

.1261 

5. 

.706 

0.113 

.0949 

.1353 

.1616 

.2000 

6 . 

.863 

0.130 

.1663 

.2093 

.2374 

.2786 




y> 

Ft 

Pressure, 

atm 

Enthalpy 
Kcal/100 gm 

Flow s 
Speed, fi 

.0854 

0.1452 

338 

2268.6 

.0827 

0.1526 

341 

1460.0 

.0809 

0.1559 

342 

1021.1 

.0783 

0.1612 

343 

653.2 

.2545 

.1381 

322 

4348.3 

.2465 

.1457 

332 

3215.0 

.2412 

.1488 

336 

2573.5 

.2336 

.1516 

340 

1823.9 

.4210 

.1367 

309 

5580.0 

.4070 

. 1403 

319 

4700.0 

.3980 

.1419 

325 

4051.2 

.3849 

. 1435 

333 

3024.1 

.5841 

,1181 

268 

8230.3 

.5638 

.1214 

293 

6700.0 

.5499 

.1226 

308 

5662 . 6 

.5296 

.1230 

324 

4151.5 

.7287 

.1163 

252 

9046.8 

.7054 

.1172 

276 

7760.0 

.6902 

.1174 

294 

6685.7 

.6680 

.1183 

319 

4663.1 

.8689 

.1194 

253 

8995.6 

.8440 

.1220 

270 

8120.1 

.8278 

.1233 

286 

7150.1 

.8040 

.1253 

317 

4895.3 



TABLE 10 (continued) 


l 

vjD 

O'* 

I 


Body 

Station 

7. 


8 . 


9. 


10 . 


11 . 


12 . 


Arc Length 
on Surface, 
Ft 

1.021 


1.346 


1.841 


2.336 


2.830 


3.325 


Standoff 

Distance, 

Ft 

x. 

Ft 

y> 

Ft 

Pressure, 

atm 

Enthalpy 

Kcal/100 

0.132 

.2352 

1.010 

.1224 

258 


.2819 

.9835 

.1254 

267 


.3123 

.9659 

.1267 

281 


.3571 

.9401 

.1283 

312 

0.177 

. 3666 

1.311 

.1268 

272 


.4253 

1.277 

.1290 

271 


.4637 

1.255 

.1300 

275 


.5200 

1.222 

.1305 

295 

0.232 

.5657 

1.7671 

.1269 

261 


.6429 

1.723 

.1292 

271 


.6933 

1.693 

.1301 

274 


.7674 

1.651 

.1306 

288 

0.300 

.7541 

2.229 

.1311 

284 


.8539 

2.172 

.1328 

280 


.9190 

2.134 

.1334 

278 


1.015 

2.079 

.1336 

284 

0.368 

.9433 

2.691 

.1260 

272 


1.065 

2.621 

.1287 

276 


1.145 

2.575 

.1297 

276 


1.262 

2.507 

.1302 

282 

0.437 

1.131 

3.154 

.1316 

286 


1.276 

3.070 

.1330 

282 


1.371 

3.016 

.1336 

280 


1.509 

2.936 

.1339 

282 


Flow 

Speed, fps 

8779.2 

8261.7 

7502.3 

5347.8 

8039.4 

8057.1 

7828.9 

6583.4 

8139.7 

8045.3 

7915.4 

7049.7 

7280.0 

7582.9 

7689.2 

7278.0 

7998.1 

7781.7 

7771.8 

7451.1 

7217.0 
7437 .3 

7549.5 

7411.8 
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TABLE 10 (continued) 


Body 

Station 

Arc Length 
on Surface, 
Ft 

Standoff 

Distance, 

Ft 

X, 

Ft 

13. 

3.820 

0.508 

1.317 

1.485 

1.595 

1.757 

14. 

4.314 

0.571 

1.509 

1.699 

1.822 

2.004 

15. 

4.809 

0.644 

1.694 

1.907 

2.047 

2.251 

16. 

5.304 

0.718 

1 885 
2.122 
2.276 
2.499 

17. 

5.798 

0.776 

2.074 

2.339 

2.507 

2.746 

18. 

6.124 

0.821 

2.194 

2.481 

2.659 

2.909 


y, 

Ft 

Pressure, 

atm 

Enthalpy 
Kcal/100 gm 

Flow 

Speed, f p ; 

3.618 

.1257 

275 

7852.5 

3.521 

.1286 

278 

7649.8 

3.457 

.1296 

278 

7653.8 

3.364 

.1303 

281 

7516.2 

4.078 

.1295 

280 

7563.3 

3.969 

.1314 

281 

7508.2 

3.897 

. 1322 

280 

7525.9 

3.793 

.1328 

282 

7405.9 

4.543 

.1282 

280 

7549.5 

4.420 

.1300 

280 

7551.4 

4.339 

.1308 

280 

7567.2 

4.221 

.1313 

282 

7433.4 

5.004 

.1265 

275 

7846.6 

4.866 

.1294 

279 

7604.6 

4.778 

.1305 

280 

7565.2 

4 . 644 

.1315 

283 

7346.8 

5 .465 

.1301 

281 

7478.7 

5.312 

.1313 

281 

7510.1 

5.216 

.1319 

280 

7525.9 

5.078 

.1327 

285 

7252.4 

5.773 

.1291 

281 

7462.9 

5.607 

.1306 

280 

7527.8 

5.504 

.1314 

280 

7545.5 

5.360 

.1327 

285 

7222.9 



TABLE 10 (continued) 


Body 
S tation 

Arc Length 
on Surface, 
Ft 

S tandof f 
Di s tance , 
Ft 

Ft 

19. 

6.281 

0.851 

2.252 

2.557 

2.740 

2.988 

20. 

6.438 

0.831 

2.311 

2.651 

2.836 

3.066 

21. 

6.595 

0.920 

2.428 

2.775 

2.951 

3.153 

22. 

6.752 

1.060 

2.628 

2.938 

3.091 

3.263 

23. 

6.909 

1.308 

2.914 

3.157 

3.273 

3.397 

24. 

7.065 

1.698 

3.356 

3.454 

3.500 

3.547 


y. 

Ft 

Pressure, 

atm 

Enthalpy 
Kcal/100 gm 

Flow 

Speed, fi 

5.920 

.1273 

280 

7557.3 

5.745 

.1298 

280 

7561.3 

5.639 

.1308 

280 

7567.2 

5.496 

.1323 

286 

7191.4 

6.068 

.1270 

279 

7608.5 

5.872 

.1293 

280 

7581.0 

5.765 

.1302 

280 

7553.4 

5.632 

.1307 

286 

7213.0 

6.326 

.1213 

270 

8106.3 

6.055 

.1219 

275 

7856.4 

5.918 

.1187 

274 

7901.7 

5.762 

.1079 

275 

7858.4 

6.720 

.1130 

256 

8863.8 

6.305 

.1062 

265 

8395.5 

6.100 

.0967 

262 

8535.2 

5.872 

.0734 

258 

8757.6 

7.169 

.0968 

225 

1030.4 

6.552 

.08613 

252 

9050.7 

6.260 

.07353 

246 

9200.3 

5.950 

.0470 

244 

9440.3 

7.686 

.0794 

184 

11954.8 

6.793 

.0657 

238 

9725.6 

6.393 

.0533 

236 

9826.0 

5.990 

.0291 

234 

9946.0 




FIGURE 19. SPHERE-CAP BODY, NON-EQUILIBRIUM SHOCK SHAPE AND 
STREAMLINE POSITIONS 
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the shock layer densities to be underestimated and the shock layer thick- 
ness to be overestimated. Since the shock layer volume is an important 
radiative heat transfer variable, the shock detachment distance has been 
adjusted in the radiation calculation to include the influence of the 
equilibrium property zone. This adjustment was accomplished by multiplying 
the existing detachment distance by a weighted average of equilibrium and 
y = 1.225 sphere-cap solutions. The error in the inviscid chemistry 
model does not significantly affect the pressure field nor the relative 
location of the streamlines. 

This error in chemistry also affects the stagnation point convective heat 
transfer and the convective heat transfer distribution. The only heat 
transfer variables. significantly influenced are the boundary layer edge 
velocity and the stagnation point velocity derivative. These variables 
have been adjusted to account for the equilibrium chemistry effects. 

The stagnation point velocity gradient is reduced by 21% from 697 sec ^ 
to 547 sec . The stagnation point pressure, on the other hand, is in- 
creased only 17 0 by the equilibrium chemistry effect. 

In Figure 20, the sphere-cap surface pressure as a function of the ratio 
of body surface distance to base radius is shown. For comparison, the 
Newtoman-Busamann and the modified Newtonian pressure distributions are 
also indicated. The Newtonian-Busemann model is a better approximation 
of the sphere-cap pressure distribution than is the modified Newtonian 
result. Both of these approximate models, however, are m error m the 
region of the corner where a local acceleration of the flow occurs. 

The flow field properties throughout the sphere-cap inviscid shock layer 
are given m further detail m Table 11. The notation follows that used 
in Table 5 and discussed in Section 4.1.1. The location of each body 
point is shown in Figure 19. Data is given along rays normal to the 
body surface. The final column in Table 11 is the enthalpy of the 
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TABLE 11 


NON- EQUILIBRIUM INVISCID SHOCK LAYER PROPERTIES ON THE SPHERE-CAP 


Body 

Station 

Arc Length 
on Surface, 
Ft 

S tandof f 
Distance , 
Ft 

i 

Ft 

| CT - 

Pressure, 

atm 

Enthalpy 
Kcal/ lOOgm 

Flow 

Speed, 

fps 

Enthalpy, with 
Radiation Loss 
Kcal/100gm 

1 .Shock 

0.04 

1.190 

-1.190 

0.0 

.161 

339 

1987.2 

339 

. 625e 



- .734 

0.0 

.153 

342 

1042.8 

325 

.375s 



- .436 

0.0 

.157 

343 

500.0 

316 

Body 



0 

0.0 

.163 

343 

110.2 

307 

2. 

0.629 

1.201 

-1.187 

0.683 

.144 

338 

2184 

338 




- .718 

0.662 

.152 

341 

1298 

324 




- .423 

0.648 

.155 

342 

807 

316 




.014 

0.628 

.160 

343 

433 

307 

3. 

1.888 

1.233 

-1.100 

2.044 

.142 

331 

3325 

331 




- .628 

1.981 

.149 

338 

2282 

320 




- .325 

1.940 

.152 ' 

340 

1771 

313 




.123 

1.880 

.155 

342 

1259 

306 

4. 

3.147 

1.243 

- .870 

3-395 

.137 

322 

4403 

322 




- .402 

3.288 

.147 

332 

3207 

318 




- .108 

3.219 

.150 

335 

2754 

309 




.341 

3.117 

.152 

338 

2105 

304 

5. 

4 .406 

1.270 

- .550 

4.700 

.133 

307 

5686 

307 




- .113 

4.583 

.138 

321 

4506 

309 




.205 

4.483 

.140 

327 

3817 

307 




.668 

4.337 

.142 

334 

2951 

302 

6 . 

5.296 

1.336 

- .288 

5.655 

.126 

294 

6690 

294 




.193 

5.470 

.129 

312 

5273 

303 




.505 

5.350 

.129 

321 

4625 

303 




.961 

5.174 

.128 

329 

3640 

302 
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TABLE 11 (continued) 


Body 

Station 

Arc Length 
on Surface, 
Ft 

Standoff 
Dis tance, 
Ft 

x, , 
Ft 

y, 

Ft 

Pressure, 

atm 

Enthalpy 

Kcal/lOOgm 

Flow 

Speed, 

f P s - 

Enthalpy , with 
Radiation Loss 
Kcal/ lOOgm 

7. 

5.665 

1.337 

- .142 

6.020 

.1225 

287 

7142 

287 




.334 

5.823 

.124 

308 

5647 

300 




.643 

5.695 

.123 

317 

4879 

299 




1.09 

5.508 

.1226 

326 

3935 

300 

8. 

5.868 

1.349 

.005 

6.348 

.1188 

275 

7870 

275 




.458 

6.297 

.1215 

303 

6040 

297 




.751 

5.935 

.1219 

314 

5116 

298 




1.176 

5.699 

.1212 

325 

4073 

300 

9. 

6.048 

1.395 

.230 

6.812 

.1153 

263 

8559 

263 




.630 

6.441 

.1165 

290 

6591 

285 




.886 

6.205 

.1183 

311 

5391 

297 




1.253 

5.865 

.1179 

323 

4289 

298 

10. 

6.163 

1.546 

.514 

7.285 

.1004 

234 

9901 

234 




.843 

6.768 

.0835 

269 

8185 

264 




1.051 

6.443 

.0659 

270 

8120 

255 




1.353 

5.982 

.0368 

244 

9466 

219 

11. 

6.300 

1.871 

.830 

7.756 

.0837 

198 

11451 

198 




1.123 

7.063 

.0671 

253 

9031 

249 




1.270 

6.635 

.0544 

260 

8677 

245 




1.475 

6.000 

.0281 

243 

9523 

218 

12. 

6,475 

2.307 

1.350 

8.288 

.068 

155 

13025 

155 




1.516 

7.370 

.0514 

234 

9916 

230 




1.572 

6.812 

.0402 

245 

9385 

231 




1.650 

6.000 

.0119 

221 

1047 

197 



TABLE 11 (continued) 


Body 

Station 

13. 


Arc Length Standoff 


on Surface, 

Dis tance, 


y> 

Ft 

Ft 

_1L_ 

Ft 

6.825 

2.850 

2.000 

8.850 



2.000 

7.638 



2.000 

6.961 



2.000 

6.000 




Flow 

Enthalpy , with 

Pressure, 

Enthalpy 

Speed, 

Radiation Loss 

atm 

Kcal/lOOgm 

fps 

Kcal/lOOgni 

.0630 

137 

13615 

137 

.0356 

220 

10566 

217 

.0249 

226 

10290 

211 

.0150 

186 

11900 

164 



fluid, corrected for radiation cooling, which is discussed in further 
detail in Section 4.2.3. The non-equilibrium species distributions are 
also discussed in Section 4.2.3. 

4.2.2 NON- EQUILIBRIUM CONVECTIVE HEAT TRANSFER 

a. Sphere-Cone , The non-equilibrium sphere-cone convective heat 

flux versus distance along the body is shown m Figure 21. For comparison, 
the equilibrium convective heat transfer distribution is also shown. The 
stagnation point convective heat transfer for the non-equilibrium case is 
slightly larger than the equilibrium result. This increase is an inviscid 
effect, arising from a higher stagnation point velocity gradient in the 
non-equilibrium case (7570 sec *), than m the equilibrium case (6340 
sec ^) . The sphere-cap vorticity interaction effect is actually smaller 
in the non- equilibrium case (1.06) than m the equilibrium case (1.12). 

The CO^ composition effect is the same for both cases. The convective 
heat transfer in the conical region of the body is slightly less for 
non-equilibrium flow than for equilibrium flow because the boundary layer 
edge velocities are slightly lower. 

b. Sphere-Cap . The non-equilibrium sphere-cap convective heat flux 
distribution is shown m Figure 22. The heat transfer distribution is 
quite flat over most of the body, but drops sharply in the corner region. 
The vorticity interaction effect for the sphere-cap case is negligible 
because of the large nose radius. 

It was found that the stagnation point velocity gradient exceeds the 
modified Newtonian result by approximately 19%; hence the sphere-cap 
stagnation point heat transfer rate is 9% larger than the modified 
Newtonian result. The effect of the CO^ composition on the convective 
heat transfer is assumed to be solely dependent on the free stream 
molecular weight and increases the convective heat transfer rate by 
approximately 8% over that for air at the same flight conditions. Radia- 
tion absorption m the boundary layer has a negligible influence on con- 
vective heat transfer. 
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HEATING RATE, BTU/FT -SEC 



SHPERE- 

CONE 

JUNCTION DISTANCE ALONG THE BODY, S/R_ 

E 


FIGURE 21. SPHERE-CONE NON-EQUILIBRIUM CHEMISTRY 
CONVECTIVE AND RADIATIVE HEAT TRANSFER 
DISTRIBUTIONS 
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HEATING RATE, BTU/FT -SEC 



FIGURE 22. SPHERE-CAP NON-EQUILIBRIUM CHEMISTRY CONVECTIVE AND 
RADIATIVE HEAT TRANSFER DISTRIBUTIONS 


-I 



4.2.3 NON- EQUILIBRIUM RADIATIVE HEAT TRANSFER 

a. Non-Equilibrium Concentrations of Significant Species . The 
non-equilibrium species profiles on three streamlines are shown in 
Figures 23, 24, and 25 as a function of distance along the streamline. 

Gas temperature and density profiles are also shown. The three stream- 
lines originate at shock angles of 89°, 79°, and 68°, respectively to the 
flow direction. The 89° and 79° streamlines correspond to the sphere-cap 
streamlines indicated by arrows on Figure 19; the 68° streamline corresponds 
to the sphere-cone streamline indicated by the arrow on Figure 17. The 

89° and 79° streamline results are also representative of the non-equilibrium 
flow in the nose region of the sphere-cone, but with a change in scale that 
correctly transforms the convective derivative (a scale change by a factor 
of about 12). 

In Figures 23, 24, and 25, the overshoot in CN concentration is seen to be 
moderate; The CN emission, Figures 26 and 27, overshoots more strongly than 
does the CN concentration because of the high non-equilibrium gas tempera- 
ture (~9000°K) near the peak in the CN concentration. The CN* concentration 
is seen to_ overshoot by a factor of about 10. The CO* concentration is not 
shown because the rates are too uncertain. The electron-density profile is 
also shown, although its contribution to the radiation heat transfer is 
negligible. However, the information may be useful m communication black- 
out studies. It is seen that a slight electron density overshoot occurs; 
the location of the peak electron density is a factor of 5 to 10 further from 
the shock than is the peak CN radiation emission. 

b. Sphere-Cone Non-Equilibrium Radiative Heat Transfer . From the re- 
sults of Figures 23, 24, and 25, it is seen that non-equilibrium gas 
properties extend over the full shock layer on the sphere-cone. 

The major radiation to the sphere-cone originates from the non-equilibrium 
overshoot of the CN red and CN violet bands. The volume emission from 
these bands along each of the streamlines noted above is shown in 
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CONCENTRATION, MOLES PER ORIGINAL MOLE 



TEMPERATURE, 10 4 °K 
DENSITY, 10~ 5 GM/CM 3 




EQUILIBRIUM 

VALUES 



STREAMLINE DISTANCE, Z, FT 


FIGURE 24. NON-EQUILIBRIUM SPECIES AND TEMPERATURE AND 
DENSITY PROFILES ALONG THE 79° SPHERE-CAP 
STREAMLINE. EQUILIBRIUM VALUES ARE ALSO SHOWN 
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TEMPERATURE, 10 4 °K 
DENSITY, 10 " 5 GM/CM 3 





i quilt bimum 
values 



FIGURE 25. NON -EQUILIBRIUM PROFILES OF SPECIES CONCENTRATION AND TEMPERATURE 
AND DENSITY ALONG THE 68° STREAMLINE INDICATED IN FIGURE 17. 
EQUILIBRIUM VALUES ARE ALSO SHOWN, 
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TEMPERATURE , 10 4 K, 
DENSITY, 'lO"-5 GM/CM^ 



VOLUME EMISSION, WATTS /CM -5 - S TER 





VOLUME EMISSION, WATTS/CM -STER 



FIGURE 27. NON-EQUILIBRIUM PROFILES OF CN RED EMISSION ALONG THE 
THREE STREAMLINES INDICATED IN FIGURES 17 AMD 19. 
EQUILIBRIUM VALUES ARE ALSO SHOWN. 



Figures 26 and 27. Behind the normal shock, the contribution from the CN 
violet band exceeds that from the CN red band by a factor of about 1.5; 
behind the 68° shock, the CN red band exceeds the CN violet band by a 
factor of about 2. The noticeable decrease in CN violet emission is due 
both to the lower gas temperature behind the oblique shock, and to an 
increase in the effect of collision limiting. 

The total radiative heat transfer to the 6 body points on the sphere-cone 
shown in Figure 17 is given m Table 12. Radiation to points on the 
rounded section between the cone and the base was not computed because 
the non-equilibrium computer program did not run satisfactorily in this 
region. 

The results listed m Table 12 include the effects of self-absorption, 
thermal boundary layer thickness, and collision limiting. Radiation 
cooling was not included because the reduction m gas enthalpy is less 
than 57 * and affects the calculated radiation by less than 27 *. The thermal 
boundary layer correction reduces the radiation by about 57„. Self- 

Zj. -f- 

absorption of CO emission reduces the radiation by as much as 6% on the 

cone surface. The magnitude of the collision limiting effect on CN violet 

emission is difficult to estimate without making calculations that omit 

collision limiting, but it is thought to reduce the total radiation to the 

cone surface by a factor of 1.5 to 2. The contribution to the total 

radiative heat transfer from the CO vibration-rotation bands and from the 
4+ 

CO band is 3 to 57„ each. 

Comparison of the non-equilibrium radiative heat transfer result on the 
sphere-cone, Table 12, with the equilibrium values shown m Table 9, 
column 3, indicates an increase by a factor of about 4 m the stagnation 
region and about 2.5 trhoughout the conical region. It is evident from 
these results that an equilibrium calculation is misleading for these 
flight conditions. 
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TABLE 12 

NON- EQUILIBRIUM SPHERE-CONE RADIATIVE HEAT TRANSFER 


Position Radiative Heat^Transf er 

(Figure 17) watts/cm 


1 

2 

3 

4 

5 

6 


20.6 

15.0 
7.2 

12.5 

22.0 

23.3 


TABLE 13 


NON- EQUILIBRIUM SPHERE- CAP RADIATIVE HEAT TRANSFER 


Position 

Radiative Heat Transfer, 

2 

watts/ cm 

(Figure 19) 

Non-Equilibrium Region 

Total 

1 

45.5 

79.8 

2 

44.6 

75.6 

3 

27.7 

54.6 

4 

13.8 

32.4 
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Comparison of the non-equilibrium radiative heat transfer on the sphere- 
cone with the convective heat transfer, Figure 21, shows that the former 
is only 107» of the latter in the stagnation region. However, on the 
conical region near the base the radiative heat transfer is 30% of the 
convective heat transfer and thus represents a significant contribution 
to the total heat transfer. 

c. Sphere-Cap Non-Equilibrium Radiative Heat Transfer . From the 
non-equilibrium species profiles of Figures 23, 24, and 25, it is seen 
that non-equilibrium gas properties cover only about 1/3 of the sphere- 
cap shock layer. The remaining 2/3 of the layer is very close to equilib- 
rium. The width of the non-equilibrium zone was found to vary very little 
throughout the shock layer, and was approximated by a constant 0.35 ft, 
shown in Figure 19. Because a major fraction of the shock layer is in 
equilibrium, the standoff distance computed using the frozen y result is 
too large by about 25%. This is corrected m the radiative heat transfer 
calculation by reducing the thickness of the equilibrium zone by 0.30 ft. 

The non- equi librium radiation emission from the CN red and violet bands 
along the .three streamlines noted above is shown in Figures 26 and 27. 

As mentioned previously the effect of radiation cooling is very small in 
the non-equilibrium region and is not included. Beyond the non-equilibrium 
region radiation cooling is included by reducing the gas enthalpy by the 
amount of the combined non-equilibrium and equilibrium radiation losses. 
This reduced enthalpy, listed in Table 11, and the gas pressure are used 
to obtain equilibrium radiation properties. Representative values of 
equilibrium volume radiation for various values of enthalpy and pressure 
are found in Table 5. 

Profiles of the volume emission with and without radiation cooling are 

o 4^ 

shown in Figure 28 for the 89 streamline. In these profiles the CO 

emission is omitted because it is so strongly self-absorbed in the sphere- 

4 •{* 

cap flow. However, the CO contribution is included in the radiative 
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VOLUME EMISSION, WATTS /CM -STER 



RIGURE 28. RADIATION PROFILES ALONG THE STAGNATION STREAMLINE OF THE SPHERE 
CAP. THE CN RED AND VIOLET BANDS AND THE CO VI BRAT I ON -ROTATION 
BAND ARE INCLUDED: THE 00^+ BAND IS OMITTED BECAUSE IT IS 

STRONGLY SELF-ABSORBED. 
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heat transfer calculation. It is seen that radiation cooling reduces the 
equilibrium volume emission by as much as 25% at the boundary layer edge. 

It reduces our calculation of the total radiative heat transfer by about 
10%, at the stagnation point to about 6% at point 4. 

The total radiative heat transfer to the 4 body points shown m Figure 19 

is given in Table 13. These results include self-absorption, radiation 

cooling, thermal boundary layer thickness. Collision limiting in the 

4+ 

non-equilibnum region is also included. The CO emission is strongly 
self-absorbed and, as a result, contributes only 4 to 5%, of the total 
radiative heat transfer. The CO vibration-rotation bands contribute about 
6 %> to the total. The remainder arises from the CN red and violet bands 
There is some self-absorption of the CN violet emission from both the 
equilibrium and the non-equilibrium regions; this reduces the computed 
radiative heat transfer by about 10%,. The thermal boundary layer correc- 
tion reduces the radiative heat transfer by about 3%. 

The contribution of the non-equilibrium region to the total radiative heat 
transfer is listed separately in Table 13 and is seen to represent 40% to 
60% of the total. If equilibrium gas properties had been assumed, the 
calculated radiative heat transfer would have been smaller by approximately 
35%. 

It is seen by comparing the results of Table 13 with the results shown in 
Figure 22 that the radiative heat transfer exceeds the convective heat 
transfer throughout the front face of the sphere-cap. This occurs because 
the shock layer is thick, increasing the radiative heating and reducing 
the convective heating in comparison to the values found on the sphere-cone. 
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SECTION 5 


CONCLUSIONS AND RECOMMENDATIONS 


Methods have been developed and described here for computing the state and 
motion variables and the convective and radiative heat transfer over the 
front face of spherically blunted conical bodies traveling at hypersonic 
speeds through planetary atmospheres under equilibrium or non-equilibrium 
conditions. The numerical results carried out for Martian entry on a 
hyperbolic trajectory show the following significant results. 

The flow field contains extensive regions of non-equilibrium gas. The non- 
equilibrium nature of the flow has very little effect on convective heat 
transfer (Figure 21), but leads to much greater radiative heat transfer 
than that computed on the basis of equilibrium flow (compare Figures 16 
and 21). 

The radiative heat transfer computed on the basis of non-equilibrium flow 
is a significant fraction of the total heat transfer, For the sphere-cap 
body shape, the radiative heating is more than twice the convective heating 
at the stagnation point (Figure 22) . Although radiative heating is con- 
siderably smaller than convective heating on the sphere-cone (Figure 21), 
this is due to the small standoff distance and the oblique conical shock 
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that occur on this shape only when the angle of attack is near zero. 

Although angle of attack calculations were not made m this study, it 
is expected that at 10° angle of attack the windward portion of the 
sphere-cone will support a shock angle and standoff distance that are 
comparable to those typical of the flow over the sphere-cap. Thus 
radiative heat transfer may exceed convective heat transfer on the sphere- 
cone at moderate angles of attack. 

The low Reynolds number associated with these calculations leads to a 
significant thickening of the boundary layer. The major result of this 
is the increase in the convective heat transfer over the value based on a 
modified Newtonian approximation. The actual convective heating to the 
sphere-cone (Figure 16) exceeds the modified Newtonian value by more than 
50% over most of the surface. This arises because the boundary layer on 
the conical portion of the sphere-cone extends beyond the low-velocity, 
high-entropy layer of fluid that passes through the spherical nose region. 

The coupling among the inviscid flow, the viscous flow, and the radiation 
transfer was found to be significant, but amenable to evaluation through 
first order corrections to the uncoupled results. 

It is recommended that m future work the non-equilibrium radiation from 
the CO fourth positive band be studied experimentally, that collision 
limiting of the CN violet band be studied experimentally, and that the 
non-equilibrium analysis be applied to the spliere-cone at small and moderate 
angles of attack in order to determine the convective and radiative heat 
transfer rates . 
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APPENDIX A 

DERIVATION OF FINITE DIFFERENCE EQUATIONS 


The left hand side and the first term on the right of Equation (6) are 

seen to be the conservation quantity averages that are desired. Let H(t) 

be the cell j+%, i+%. The area of the cell is denoted at the start of the 

time step by G , , , and is determined from 

J "H , i+5 


G. , . , = ff dydx 
o 

Similarly, the area of 0(t +T) is denoted by 


(Al) 


dydx 

n(t +t) 


(A2) 


The mean value of the conservation quantity f in the cell j+^, ±+k at 
t +T is defined by 


0(t +T) 

= o 


Jj f(x,y,t +T)dydx 


j+%,i+% 


(A3) 


with a similar definition at t for f . ,, ... 
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The boundary motion gives a contribution only on the segments o£ the 
boundary j, i+% and j+1, i+k. The ray portions of the boundary are 
stationary and have = 0. The boundary motion contribution to the 
change in f is, therefore. 


t +T 
o 

J* dt If f(x,y,t)V ds 

t o o0(t) 


- f . . ,1 G . j 


^ j +1 , i+k^ j +1 » i+^ 


(A4) 


G^. ^ +1 denotes the volume swept by the boundary, j, i-f% as it moves inward. 


and 


t +*T 
o 


G j ,i+% 




(j ,i+%) 


V ds 

n 


G 


j+l,i+i> 


t +T 
o 

= J dt 

t 

o 


f V ds 

(j+l,i+%) n 


(A5) 


(A6) 


This definition corresponds to taking the uniform convention for the 
boundary velocity throughout the grid as positive when it moves toward 
the body. The velocity is taken outward from each cell, and is in- 
convenient for use in moving from cell to cell as it must be continually 
redefined. 


id-*? i+% i id% 

Since Equation (A4) is already of order T compared to f G J ’ 

the value of f. need only be correct to zero order. 


f . . - %(f. , o + f.o, IjA ) + 0([f. . .]) 

J,lH J-%,1+% J+%,l+*2 


where the bracket denotes the jump in f at the cell boundary. In Godunov's 
method a specific form is assigned to the first order term in this expression 
for all boundary quantities. The resulting first order terms stabilize 
the calculation. 
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The contributions from the line integral of P,Q about the cell boundary 
are seen to have the -form t(f,S) for each edge. Here F is taken to be 
the two component vector IP,Q] and S the two component vector with 
magnitude equal to the length of the edge and direction along the out- 
ward normal to the edge. The notation (F,S) is the scalar product 
between the two vectors. S can be evaluated for the edge at time 
for first order accuracy in Equation (6), since the difference of T(F,S) 
for opposite sides of the cell is already of first order relative to 
£j + 2 3 + i-H^ •j'hg fi ow quantities appearing in (F,S) must be 

evaluated by Godunov's method for stability. Finally then, the flux 
contribution is seen to be 


t +T 
O 

J* dt f Pdy - Qdx = T I (F,S) 
t oQ(t) (edges) 


(A7) 


For axisymmetric fiow the source term in Equation (6) remains to be 
approximated. Again it is of first order compared to the leading term 
f J 2} 2 G j 25 2 so that the geometric factors and flow variables can be 

evaluated at t = t . When y {= 0 inside the cell, this term will be 
approximated by 


where 


and 


t +T 

o 

J dt J7 “ R(x,y,t)dsdy 
t o H(t) y 


TR . . i - 1 B , ,j i 


- £i' i 




R 


j +h,±+h i+%^ 


(A8) 


(A9) 


(A10) 


For cells bordering the axis, the expression (A9) is not convergent and 
must be replaced by the convergent integral 
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-If ?<>,<•* 

n(t o ) y 


and a compatible modification of R to R 


o 


R°(f) = 


p 

Pu 

Pv 

p(e +p 


(All) 


(A12) 


Denote Godunov's evaluation of p, u, v, p, e on the cell boundaries by 

R, U, V, P, E respectively. Let W denote the fluid velocity normal 

to the segment j,i+k directed toward the body and computed from V,U. 

Similarly let W. , be the fluid velocity computed from V.U normal to 
3+%>i 

j+^,i and directed along the shock layer from the axis. By Equations 
(A3), (A4) , (A7) and (A8) the finite difference approximations to 
Equation (6) are 


Conservation of Mass 

VM-i - T[ » ws > J+ i,ui - (raB) j*.i ] 


-LE(TWS«), +1j1+% - R(TWS-H5) Jjl+i ] - Tk(pvB) j+t _ i+% ) 


(A13) 


Conservation of x momentum 

» J+ * 14 * - [(PGP^^VUPUG) . . 

-T[(EUWS) j+iii+l - 

-[R0<THS«) j+1j144 - RUOHS+G). ^] 


-T[(PSsil,o) j+1>i+4 - (PSsxnc)^.^] 
-T[(PS 5l „9) j+%il+1 - (PS.x„8) j44il ] 


(A14) 
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Conservation of y momentum 

vi+i.l+i , C(p(3) ^,^]-l t ( p vG) . 

J+'gji+'S 


-t[rvws) 


J+%,1+1 




~[RV(TWS+G) • , 1 , - RV (TWS+G) 'j 

j +1,1+3 

-t[ (PS cosB) , - (PScosQ) , j 

J+%>1+1 

+r[(PScosa) j+lii+% - (PScosc) ji+% ] 

-Tk(pv 2 B) j+%ji+i } 

Conservation of energy 

e 3+%> i +% , c (PG ,^. i+ 4j-H(p E G) . i+t 

j+3,1+3 

-T[(EEHS ) J+4>1+1 - (REHS) j+5s i ] 
-CrE(tWS+G) . , - RE (tWS+G) . .] 

j+l,i+% j 3 i+% 

- tC( ™ s) j H >1 h - (PWS) J+ i, 1 ] 

-T[(PWS). + i ;i+i - (PWS) j>iH ] 


-Tk[pev+pv)B] . J 
J +g , 1+3 


(A15) 


(A 16 ) 

The difference Equations (A 13 ) through (A 16 ) contain the geometric quantities 


^i+h, i+% . 


5 G j+^,i +V G j 3 iiV S J 3 r+%’ S j+%,i 5 B j-^,i+% lf y ^ ° and 


B 


j +%,!+% 


! when y = 0 and sine . , , cosa 




J , 1+3 


Thus the area of the cell is explicitly 
G. 


3. ,i ,i “ fl dydx = 0 i - 0. , . j 

J+3,i+3 y J J,i+ J 2 J+1,1+% 


(A 17 ) 


where 

0 


- 2 2 

= ^ cos9 .^sinQ _ - £ cosQ,sin9 + Oz. ~ -t, .cos9 f ,, 

’ " N ^l+1 J , 1+1 1+1 


= Ht 

J,i+% 2 J > i + l i+l 1+1 J ,1 1 1 


-§ . + £ .cosQ )(£. ^sinQ. ,, + t .sin0 )] 

i J,i i J,i+1 i+l J,i i 


S lmi larly 


G J+%,i+^ = JJ dydx = 0 j,i+% _ 0 j+-X 3 i+% 


(A 18 ) 

(A 19 ) 


fi(t 0 +T) 
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with 

$j> 1+ ^ = %[ ’ 1+1 ) 2 cos0 sm9 - (^ J,1 ) 2 cos9.sin6 

i+l l+l 11 

+ (5 , - 5 1 ^"^cos9 . . - £. + -L~* ’ ^cosG . ) (-t~® 5 1+ ^sin9 , + 5 ^sin9 . ) ] (A20) 

' 5 1+1 X+l 1 ' 1+1 

where 5 1 is the location of t .at the end of the time interval. The 

3 

area swept by the boundary j,i+^ is 


G j ’ 1+% = 0 3,i+h _ 0 


J»i+% 


(A21) 


With the same definition holding when j is replaced by j+1. The lengths 

of the sides S“ S , are given by 

J, i+% J+fc,i * 


S. , = l. . - l. , . 

j+%,i J,i J+l,i 


(A22) 


and 


S . ,i = C(§ ,, - l cos9 . - § + t . cos9 ) + (l sinO 

J,i+is i+l j,i+l i+l i j,i i J,i+1 i+l 


2i% 


-■L . sm9.) ] 

J = i i 


(A23) 


The sme and cosine of the slope of a ray are determined from 9^ while 
the sine and cosine of the slope of a segment on the transverse inter- 
polatry polygon are 


t sm9 , , - -C . sm9 

l,i+l i+l 1,1 l 

sino. , = — ■ u Li 

J > 1+2 S , I 

J , l+% 


(A24) 


and 


cosa, . , 
J>i+2 


5 n-i - iagiH ~ ^i.i 0056 , 

S J j i+^ 


(A25) 


The factor of the source term for axial symmetry, when y 0, is 
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*j+%,i+% n(t y ) dydX I j +} 2 > :L+1 + I j+l,i+ J 2 


I.,» . - I. 

J+%,1 J,l+^ 


(A26) 


where 


J+,>1 1 Vi,i 


(A27) 


I. a = t§ ,, - •£. (cos© - sin0 tana. , )] 

j , \+% l+l j,i+l i+l i+l j,i+% 


...sinG 
.1 ,!+!_ - i+l N 

^ ( t .sinG. ' ) 
J»i i 


(A28) 


If y = 0, the cells border the axis of symmetry, and B i i is given by 

J+^>"2 


B ° + ^ % = J*/ ^ dydx 

J 2,2 d(t y 

o 


(A29) 


Since v = 0(y) as y “* 0 the above integral is convergent. The integrand 

in Equation (A29) can be replaced by 3v/by| to first order accuracy 

y=0 

throughout the cell. An integration of Equation (A29) can then be performed 
to yield 

Sv , 

(A30) 


Vm ' f d * sT <y u (x) " y -t (x)) 


where y (x) an d y^(x) are the upper and lower boundaries of the cell. 

All cell boundaries are straight lines so that y and y, are proportional 
to x. Moreover, let 3v^/oy be approximated by the mean value of V on 
the boundary segment as determined by the Godunov scheme, divided by the 
mean y coordinate, uy, of the segment. Therefore 


and 


ov V . , 

= -1: + 

^2>y Ay ^ 

J 


(A31) 
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(A32) 


dv 

( — -) 

j+%,i 

From Equation (A30) it is seen that 




1 

1 


V j 3 % AX i,^ + V j+%,l Ax j+% 3 l V j+l,i; AX j+l > 


with the definitions 


(A33) 


Ax. i = 5 ■, - . cosQ. - 5 + •£■ 

3,h 1 J,1 1 ° 3.° 


Ax j i = (^- i " -^-.7 ,)cos9 
J+%,1 3,1 j+1,1 l 


(A34) 

(A35) 


Equations (A33-A35) along with the modified definition of R m Equation 
(A12) permit the product RB to be evaluated for cells adjoining the axis 
of symmetry. The definitions of Equations (A26-A28) along with R given 
in Equation (4) permit evaluation of the product in all other cells. 
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